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Over the last 20 years, supernovae have become a key tool to constrain the expansion history of the Universe 
through the construction of Hubble diagrams, using luminosity distances to supernovae belonging to the "la" 
subtype. This technique was key for the discovery that the expansion of the Universe is now accelerating. We 
review the principle and difficulties of the measurements, the classification and diversity of supernovae, and the 
physics of the explosion. We discuss the systematic uncertainties affecting the cosmological conclusions with 
some emphasis on photometric calibration. We describe the major supernova cosmology surveys, the presented 
analyses and their conclusions, together with the present status of the field. We conclude on the expectations for 
the near future. 
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1. INTRODUCTION 

Kinematics of test masses probe the sources of gravitational 
field, both in a Newtonian framework and in general relativ- 
ity (GR). In a Universe following the cosmological principle, 
namely homogeneous and isotropic, the motion of test masses 
(e.g. galaxies) constrains the mean density of the Universe 
and its evolution with time. 

In 1927, Lemaitre derived from Einstein's general relativity 
equations that the velocity of recession is proportional to dis- 
tance to the first order, and also extracted from the measure- 
ments of distances to galaxies (deduced from their apparent 
brightness) and their recession velocity (from their spectral 
shift) that this property might be realized in nature [l] 1 . In 
1929, Hubble presented distances to a sample of galaxies ex- 
tending further out and also proposed a linear relation between 
distance and recession velocity [3]. This was the discovery 
of the expansion of the Universe and the Hubble law states 
that distance is proportional to recession velocity, at least at 
a small velocity, and ignoring "peculiar motions" due to spa- 
tial density variations. The proportionality coefficient is now 
called the "Hubble constant" Hq. Although its value has been 
revised by almost an order of magnitude since Hubble's ini- 
tial determination, the discovery of the proportionality itself 
remains one of the major early steps of observational cosmol- 
ogy. In modern parlance, the fractional spectral shift towards 
the red, the "redshift" z has replaced the recession velocity: 
Z = /Ireception/ ^emission ~ 1- Because it is a direct observable, 
redshift is used by cosmologists to index cosmic time. In the 
current standard model of cosmology, z — 0.1 corresponds 
to ~10% of the current age of the Universe (about 13.7 Gy), 
z = 0.5 to 40%, z = 1 to 55%. The first stars were probably 
formed around z = 10 (0.5 Gy after the Big Bang). 

The Hubble constant Hq has the dimension of the inverse 
of a time, usually expressed in units of km/s/Mpc. It defines 
the expansion rate today, and sets the time and distance scales 
of the expansion. Without diminishing the importance of the 
Hubble constant determination, one should note that its value 
does not convey sufficient information to define the current 
average energy content of our Universe: GR relates the ex- 
pansion rate at a given epoch with energy density and spatial 
curvature at the same epoch. The Hubble constant defines the 
average density of the Universe only if the latter is flat (i.e. has 
zero curvature). The density for which the curvature changes 
sign is called the critical density and depends only on the Hub- 
ble constant and Newton's constant. 

In order to constrain gravitation sources driving the expan- 
sion, one has to measure the expansion beyond first order in 



Note that the English translation [2], probably carried out by Lemaitre him- 
self in 1931, misses the key paragraphs that propose the observational evi- 
dence for expansion. 



recession velocity (or redshift). Namely, as one studies the 
velocity variations with time of free falling bodies to measure 
the gravitational field of the Earth, one should probe how the 
expansion rate evolves with time in order to characterize the 
mean Universe content. In the GR framework, the expansion 
rate evolution provides a handle on the fractions of various 
Universe components, including curvature. Practically, prob- 
ing the distance-velocity relation beyond first order makes it 
possible to constrain the evolution of the energy density in the 
Universe, which in turn constrains the proportions of various 
components in the admixture, after assuming how each of the 
involved components dilutes with expansion. We will discuss 
in §2.3 why determining the proportions of the various com- 
ponents does not require measuring the Hubble constant Hq. 

How practical is it to measure the distance-redshift rela- 
tion? Measuring a redshift might be considered as relatively 
easy, at least if the target is bright enough to allow one to ac- 
quire a spectrum: redshifts of galaxies are routinely measured 
to three decimal places from spectroscopy in the visible. But 
measuring cosmological distances is not as straightforward. 
Hubble derived distances to galaxies from assumptions about 
intrinsic luminosities of stars. The overall intrinsic luminosity 
scale turned out to be wrong and yielded a Hubble constant 
about 7 times larger when compared with modern determi- 
nations. But the real weakness of galaxies as distance indi- 
cators is in fact more serious (e.g. [4]): these are evolving 
too rapidly with time to provide reliable distances, even if one 
concentrates on relative distances. Supernovae were then in- 
troduced and used as a more reliable distance indicator [5, 6], 
and eventually proposed to constrain the evolution of expan- 
sion rate [7] . 

Walter Baade and Fritz Zwicky introduced, in 1934, the 
term "supernova" to describe extremely bright explosions of 
stars that appear and fade on timescales of months to years. 
Supernovae are classified today into a handful of types and 
subtypes. The subtype of type la supernovae (SNe la) displays 
reproducible luminosities: at low redshift, the apparent bright- 
ness scales as the inverse square of the redshift (e.g. [5, 8]). 
By extension, plots of apparent brightness of similar objects 
versus redshift are often called Hubble diagrams. 

Precision supernova photometry started in the 1990s with 
the advent of CCD imagers mounted on telescopes. The 
Calan-Tololo survey delivered, in 1996, precise and homoge- 
neous CCD-based photometry in several bands for 29 SNe la 
events [9], and derived distances to those [8]. Although the 
Calan-Tololo survey should be regarded as the keystone of 
supernova cosmology, redshifts were limited to z < 0.1, too 
low to constrain the distance-redshift relation to second or- 
der. At about the same time, two teams targeted the detection 
of SNe la at significantly higher redshifts, also using CCD 
imagers. After a slow initial period, both teams eventually 
succeeded in finding and measuring redshifts and distances 
to about 50 SNe la events in total, at (z) ~ 0.5. Both teams 
compared their own distant sample to the Calan-Tololo nearby 
events (§11.2 & [10, 11]) and reached, in 1998, the same 
striking conclusion: the distant events appear too faint com- 
pared with nearby ones, in any matter-dominated Universe. In 
terms of kinematics, this means that the expansion is accel- 
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erating at late time, which is impossible under matter domi- 
nation (where matter is defined as "tracing the expansion"). 
Cosmologists named the conjectured fluid that sources the ac- 
celeration "dark energy", and it represents about 3/4 of the 
total energy density now, so that our Universe is no longer 
matter-dominated. The discovery of accelerated expansion (or 
dark energy) constitutes one of the major recent discoveries in 
physics: it causes persistent perplexity (for reasons summa- 
rized in §2.5), and still sources an enormous flow of scien- 
tific literature. The discoverers of accelerated expansion were 
awarded the 2011 Nobel Price in Physics. Dark energy stud- 
ies constitute the main scientific driver for large astronomical 
projects such as LSST and Euclid (§12). 

There had been earlier hints that our Universe might not 
follow the flat matter-dominated scenario. Different indica- 
tions of a low matter density (i.e. well below the theoretically 
favoured critical density) were presented in the early 1990s, 
mostly based on galaxy clustering (e.g. [12]) or primordial 
nucleosynthesis and galaxy clusters (e.g. [13]). The argu- 
ments in favour of a low matter density turned out to be cor- 
rect, but this does not directly imply an accelerated expansion. 
In the 1990s, settling between a (observationally favoured) 
low-density Universe and a (theoretically favoured) critical 
Universe constituted a key question, and the kinematics of ex- 
pansion seemed an obvious decisive observational test, if pos- 
sible. The enterprises aiming at producing Hubble diagrams 
of SNe la were conceived within this framework. 

Within the cosmological principle and GR, acceleration of 
the expansion happens when a fluid that does not expand (or 
expands slowly) starts to dominate the energy density budget. 
The simplest cosmological model that describes the observa- 
tions consists in a flat Universe consisting, today, of about 
1/4 of matter (diluting with expansion and time) and 3/4 of 
"dark energy", insensitive to expansion and of constant den- 
sity. Their densities were equal at z ~ 0.45, and matter domi- 
nated earlier, while the initial domination of radiation energy 
over matter lasted until z - 3300. A possible incarnation of 
dark energy of static density is the cosmological constant A 
introduced by Einstein in the GR field equations. It plays, in 
the cosmological framework, exactly the same role as a uni- 
form fluid with static density. The "ACDM" cosmological 
model consists (at late time) in this admixture of a fluid of 
static density A and matter which has to be dark and cold, 
hence coined cold dark matter or CDM. 

Cosmologists have considered a broader range of dark en- 
ergy models and describe them in terms of the "equation of 
state" of dark energy, namely how the dark energy density 
varies with time and redshift (§2.1). Describing dark energy 
using the equation of state concept should not be taken too lit- 
erally: this constitutes a phenomenological description which 
does not require the assumption that acceleration is indeed due 
to some fluid. This description of the phenomenon applies to 
actual fluids, to the cosmological constant, and it can even be 
rephrased in terms of alterations to GR if one ignores density 
perturbations. 

The discoverers of the accelerated expansion from super- 
novae readily considered alternatives to the cosmological con- 
stant paradigm, and found their observations to be compatible 



with it(§ 1 1 .2). Since then, the result has been consolidated 
mainly on two fronts: firstly, second-generation supernova 
surveys have been conducted and have confirmed the early 
results with more events and better measured distances ex- 
tending to higher redshifts. Secondly, we now have strong 
evidence for dark energy from other cosmological probes al- 
though with less sensitivity (e.g. §4.1 in [14]). To date, large- 
scale cosmological observations have been described well by 
the ACDM model (e.g. [15]), sometimes called the concor- 
dance model. 

Constraining cosmology from distances to supernovae re- 
lies on a strong hypothesis: SNe la are good "standard can- 
dles". The tight scatter of the Hubble diagram of SNe la in- 
dicates that the events are reproducible over a small redshift 
interval, but this does not prove that distant and nearby events 
are identical. They could, indeed, be of different nature be- 
cause distant events happen earlier in star formation history 
than nearby ones. This whole issue is referred to as "evo- 
lution" and can be addressed observationally: the standard 
candle hypothesis can be falsified, and is reinforced by ev- 
ery failed attempt to do so. We will, however, see that the 
picture we have of the SN la explosion mechanism makes it 
a fundamental "self-calibrated" process, weakly sensitive to 
environmental influences. 

Photometric calibration has recently become a serious issue 
for supernova cosmology. Although the core of the method 
consists in comparing fluxes, these fluxes have to be mea- 
sured in different observer bands because of redshift: nearby 
supernovae are usually measured in blue and green bands, and 
the same restframe region lies around 1 fim for supernovae at 
z ~ 1 . To take full advantage of current supernova statistics, 
these flux ratios should be calibrated to better than 1%, which 
hits fundamental limitations of current astronomical photo- 
metric calibration techniques. Today, supernova cosmology 
is probably the most demanding field of astronomy in terms 
of photometric calibration across the visible and near infrared 
range. 

The plan of this review goes as follows: §2 covers the min- 
imal necessary cosmological background, how distances con- 
strain cosmological parameters, and a summary of the leading 
hypotheses for dark energy; §3 introduces the minimal back- 
ground in astronomy techniques required for this review, for 
readers unfamiliar with the field; §4 covers the observational 
aspects of supernova taxonomy; §5 concentrates on the "la" 
subclass delivering almost all cosmological constraints, and 
§6 covers the pre-explosion and explosion issues. We then 
enter the observational and reduction techniques (§7), and de- 
vote the next section (§8) to light curve fitting, almost always 
involved in distance estimation. §9 covers photometric cal- 
ibration, the basement of supernova cosmology, and studies 
of other astrophysical systematic uncertainties are reviewed 
in §10, including gravitational lensing (§10.2). §11 covers 
supernova surveys and their cosmological conclusions; §12 
briefly discusses the prospects of the field. We summarize 
and conclude in § 13. 

For interested readers, we provide here an incomplete set 
of reviews with similar or highly related scopes: [16] parallels 
this one; [17] concentrates on explosion models; [18] focuses 
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on the physics of SNe la in relation to their cosmological use; 
[19] discusses the determination of the Hubble constant which 
we will not discuss here; [20] presents theoretical and obser- 
vational aspects of dark energy and the accelerating Universe, 
while [21] studies dark energy probes in great detail. 



2. COSMOLOGICAL FRAMEWORK 

2.1. The cosmological principle and Friedman equation(s) 

The cosmological principle states that the Universe is 
isotropic and homogeneous. Although we have ample evi- 
dence that this is wrong on small distance scales, the quasi- 
isotropy of the cosmic microwave background (CMB) [22], 
and the homogeneity of the distribution of galaxies [23] sup- 
port this framework for large scales. Within General Relativ- 
ity, this principle is encoded into symmetries of the metric: 



ds 2 = dt 2 



■ R 2 (t) 



dr 2 
1 -kr 2 



+ r 2 (d0 2 + sin 2 6d<p 2 ) 



usually called the Friedman-Lemaitre-Robertson- Walker 
(FLRW) metric. R(t) is called the scale factor, and k = -1, 
or 1, is the sign of spatial curvature. We set "c = 1" in 
what follows. Objects with constant coordinates (r, 6, <f>) are 
called comoving. Cosmologists attach comoving coordinates 
to the matter fluid, so that R(t) scales with the separation of 
matter tracers (typically galaxies). In the FLRW framework, 
it is easy to show that photons (following paths with ds 2 = 0), 
emitted by comoving sources and detected by comoving ob- 
servers, see their wavelength scale with R(t): 
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where z is called the redshift of the emitting source. GR postu- 
lates relations between the metric and the gravitation sources. 
For the FLRW metric, these are two differential equations for 
R(t), the Friedman equations [24]: 



h(z) = [rJ =— p -RHT) 



R 
R 



-^(P + 3p) + f (2) 



where A is the cosmological constant, p stands for the proper 
energy density and p for the pressure. H(z) is the expansion 
rate and Hq = H(z = 0). The energy conservation equation: 



dt 



{pR s ) = -3pR z R 



(3) 



relates pressure to density evolution and also applies sepa- 
rately to the various fluids in the Universe. For non-relativistic 
matter, pmR 3 is constant and hence pu — 0. A fluid with 
static density (p — 0) has p = -p. The cosmological constant 
A plays exactly the same role as a fluid with static density 
Pa = — pa = A/8^G, and can be summed into the p and p 



quantities in the Friedman equations. Relation (3) can be ob- 
tained by eliminating R between equations 1 and 2. 

Fluids can be characterized by a relation between p and p. 
The equation of state of each fluid wx is defined by px = 
wxpx, and equation 3 implies p x (t) oc R(t)~ 3 ^ 1+Wx \ Given the 
densities at one epoch (e.g. now) and the equations of state 
of the fluids of the Universe, one can solve the first Friedman 
equation (equation 1) for R(t). The critical density is defined 
as the energy density for which k — 0: 



Pc = 



3_Hl 
8nG 



where Hq = (R/R)„ ow is the Hubble constant. Densities are 
usually parametrized through their value "today", in units of 
the current critical density: 



D.M = — - 



Pm 8^Gp A 



3H 2 



Qa = 



3H 2 ' 



Clt = - 



and the first Friedman equation implies 1 = Q.m + Qa + ^a- 
This assumes that our Universe is composed only of matter 
and A. One might consider other fluids such as an implemen- 
tation of dark energy different from the cosmological constant 
(see §2.5) and radiation density. 



2.2. Cosmological distances and dark energy 

Cosmological distances should be defined from observ- 
ables. For example, one defines the luminosity distance 
through the energy flux / of a source of proper energy lu- 
minosity £. through d 2 L = £.IAnf. For a source observed at a 



redshift z, it reads [25]: 



d L (z) 
H 2 (z) 



(1+z) 



Sin { |Q{ 



|l/2 



r d Z ' 

Jo H(z')/H Q 



(4) 



[Q M (1 + z'? + O a + n k (l + z'f + Q*(l + z) A \ (5) 



where S in(x) = sin(x), x, sinh(x) for k = 1,0,-1; note that the 
expression is continuous in = 0. In z = 0, the integrand is 
1 and the Taylor expansion reads diiz) — z/Hq + 0(z 2 ), so that 
densities are not constrained by distances at a small redshift. 
Distances to higher redshifts become sensitive to densities, al- 
though in a somehow degenerate way. For sake of complete- 
ness, the expression for H 2 (z) (5) contains a radiation term 
Q.r, which now represents less that 10~ 4 of the energy budget 
and is usually neglected. It represents about 10% of the total 
atz ~ 1000. 

One can generalize equation 4 to alternatives to A by re- 
placing the £2 A term by the (reduced) density of the consid- 
ered fluid. For a constant equation of state, Q\ — > £lx{l + 
z) 3< - \ and for a varying equation of state w(z), Qa 
Q x exp[3f zl -^dz']. 

The discovery of accelerated expansion relied on super- 
novae at z — 0.5 being observed to be fainter (i.e. at a larger 
distance) than expected from their low redshift counterparts 
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in a matter-dominated Universe. Inspecting equation 4, one 
notes that for a flat Universe {Q.k = 0, Q.a — 1 ~ QmX a posi- 
tive £2a indeed increases distances (at given z) with respect to 
Qa = 0. The argument also holds for non-flat universes. 

The angular distance iJa is defined via the apparent angular 
size 6 of an object of comoving physical size D : <Ia = D/9. 
Because photons follow null-geodesies of the metric [26], we 
have di = (1 +z) 2 aU and hence di and d^ convey the same cos- 
mological information, di measurements rely on "standard 
candles" while d& measurements rely on "standard rulers". 



2.3. Standard candles and Hubble diagram 

Standard candles are defined as objects of reproducible lu- 
minosity and enable one to measure distances through their 
observed energy flux: / = £,IAnd 2 L {z), where X. stands for 
the proper intrinsic luminosity (energy per unit source proper 
time). Practically, we have to measure fluxes within bounded 
spectral bands, and in order to compare objects at different 
redshifts we will do our best to compare the same restframe 
bands which implies measuring high-redshift objects through 
redder (observer) bands than nearby objects. In supernova 
parlance, a Hubble diagram reports the apparent flux (or mag- 
nitude) of standard candles as a function of redshift. 

If one measures a Hubble diagram, we note that all fluxes 
scale with the factor £Mq, which conveys no information 
about cosmological densities. In the cosmological analysis, 
this overall factor will be estimated together with the other 
cosmological parameters: it enters as a nuisance parameter 
in the analysis. Since Hodiiz) does not depend on Ho, the 
Hubble constant enters the problem only through this global 
flux scale. If the luminosity of the standard candle _£ can be 
estimated without using the flux-redshift relation, measuring 
-CHq from the latter constrains Hq. We hence note that ob- 
jects with an unknown luminosity do not help determining Hq 
from a redshift-flux relation, and the initial determinations of 
the Hubble constant were off because a wrong luminosity was 
assumed. We also note that determining Ho entirely decou- 
ples from measuring the reduced densities from a Hubble di- 
agram, and we refer to [19, 27] (and references therein) for 
a detailed description of the role of SNe la to measure Ho- 
We just mention here that, using the standard candle property, 
SNe la are used to extrapolate distances from local events to 
distant ones, and constitute the last step of the "cosmic dis- 
tance ladder" used to establish the value of Hq and hence the 
cosmic distance scale. 



2.4. What can we probe with distance data? 

A set of distances as a function of redshift probes a 
weighted sum of densities (equation 4), with weights evolv- 
ing with redshift. Since distance measurements deliver a sin- 
gle number at every redshift, a redshift span is mandatory to 
break degeneracies of cosmological parameters [28]. As dis- 
cussed in the above section, supernovae impose (at least) one 
extra nuisance parameter, which illustrates why measuring a 



single cosmological parameter like Q.m in a matter-dominated 
Universe requires a finite redshift interval. 

Although distances over a finite redshift interval formally 
break the degeneracy of the sum, supernova surveys still de- 
liver elongated constraints in the (Qm, Qa.) or (Am> w) planes 
(see e.g. [10, 11]), even for surveys with the largest redshift 
span (e.g. figures 8 and 9 in [29]). 

Accelerated expansion is indeed supported by supernovae 
alone, since the first measurements (discussed in § 1 1 .2) and 
even more today (see figure 8), but with some caveats, (e.g. 
[30]). However, distances to supernovae are most usually 
merged with other cosmological probes in order to produce 
the tightest cosmological constraints and to decouple param- 
eters. Standard examples include anisotropies of the cosmic 
microwave background (see e.g. [15, 31]), and baryon acous- 
tic oscillations (BAOs, see e.g. [32, 33]) measured in the cor- 
relation function of galaxies. BAO mostly deliver an (angu- 
lar) distance constraint, with the same degeneracies as lumi- 
nosity distances, but the size of the measured "standard rod" 
also depends on cosmological parameters (mainly Q.m in the 
concordance model, see e.g. [34]), and the constraints hence 
exhibit different degeneracies than mere distances. As shown 
in [33], BAO now deliver distance-redshift relations similar 
to the ones obtained from SNe, but not yet as precise. CMB 
anisotropies probe the same geometrical feature as BAO, but 
at a much higher redshift (z - 1089) which makes CMB con- 
straints even more complementary to distances to supernovae 
typically at z < 1 to 2. 

Even if supernovae on their own constrain only combina- 
tions of cosmological parameters, they are currently key to 
constrain dark energy: ignoring supernova distances degrades 
dark energy constraints by factors of ~ 2 or more (see e.g. 
table 4 and fig 12 of [15]). 

2.5. A few possibilities for dark energy 

Right after the discovery of the accelerated expansion, 
modellers produced an enormous corpus of literature devoted 
to implementations of dark energy. Those can be classified 
into three categories: the mere cosmological constant, vac- 
uum energy, and alterations to GR. All these implementations 
involve new physics. A fourth class of proposals for dark en- 
ergy, free of new physics, consists in noting that the effective 
Friedman equation in a non-homogeneous matter-dominated 
Universe plausibly contains a term somehow mimicking dark 
energy, increasing with density fluctuations. We now briefly 
describe these possibilities and refer to [35] for more details. 
From an observational point of view, the current tightest con- 
straints (§11, 11.5) are compatible with a static dark energy 
density, which can be incarnated by the cosmological con- 
stant. 



2.5.1. The cosmological constant 

In the original papers that introduce GR, the cosmologi- 
cal constant was absent. It was later introduced by Einstein 
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himself in order to allow for a static Universe (see [36], and 
[37] for the historical aspects). Although the addition became 
useless after the discovery of expansion, it is very natural to 
allow for such a term in the Einstein equations. One even ex- 
pects such a term because it is algebraically indistinguishable 
from the contributions as gravitational sources of all "zero 
point energies" (also called vacuum energy) expected in the- 
ories of particle physics. "The cosmological constant prob- 
lem" (see [37] and references therein), is that the expected 
energy density from zero point energies of particle physics 
theories is "about" 120 orders of magnitude larger than what 
astronomical observations indicate 2 . Several classes of solu- 
tions together with their weaknesses are discussed in [37]. For 
example, super-symmetry would solve the problem because 
fermions and bosons then contribute equally but with opposite 
signs to vacuum energy, but exact super-symmetry also can- 
cels the anomalous magnetic moment of the electron, which 
is now measured to be non-zero to nine decimal places [38]. 

So in practice, one is forced to assume that the contribution 
of vacuum energy to gravitation sources is zero or extremely 
small, for some yet unknown reason. This does not forbid the 
introduction of a fundamental cosmological constant directly 
in the theory, and the observational consequence is that the 
associated (dark) energy density does not vary with time: it 
has an equation of state w = -1. In addition to postulating 
one more fundamental constant, this incarnation of dark en- 
ergy raises the "coincidence problem" (e.g. [35]): how come 
that the matter density varying as (1 + z) 3 and dark energy of 
constant density now take similar values? 

2.5.2. Quintessence 

The word "quintessence" was introduced in [39] to name 
a (fifth) component (after baryons, dark matter, neutrinos and 
radiation) in the Universe, proposed there as a scalar field, 
with a self-interaction parametrized with its potential energy. 
In most implementations, it interacts with other fluids only 
through gravitation. The dynamics of these models depend on 
the initial conditions and the shape of the potential, for which 
we have little theoretical guidance. In order to be less arbi- 
trary than the cosmological constant paradigm, these models 
try to avoid the "coincidence problem" by proposing a po- 
tential in which the scalar field density tracks the matter den- 
sity, being subdominant at early times and becoming slightly 
dominant at late times, and where the dynamics exhibit an 
attractor, making the system insensitive to initial conditions. 
Power-law and exponential potentials (e.g. [39^43]) deliver 
dynamics roughly compatible with the observations. Such po- 
tentials not only emerge in effective field theories (e.g. from 
extra dimensions or standard model extensions like SUGRA), 



but also provide the small numerical factors required to gener- 
ate the very low observed density [44]. The equation of state 
of a scalar field is confined to [-1,1]: at w = — 1 the poten- 
tial energy dominates, while the kinetic energy dominates at 
w = +1. In order to reproduce the observational constraints 
(i.e. w ^ -1), quintessence dark energy models should be es- 
sentially static at late times. "Natural" quintessence models 
commonly exhibit late equations of state around w >~ -0.8 
(e.g. [44]) which falls about 3 cr away from the latest mea- 
surements. Without fine tuning, the effective mass of the 
quintessence field is of the order of Hq [45] (i.e. 10~ 33 eV), and 
hence its spatial inhomogeneities only develop on scales com- 
parable to the horizon today. Well-motivated quintessence 
models should "protect" such a small mass from quantum cor- 
rections by some symmetry (as proposed in e.g. [43]). 

As measurements increasingly favour a static dark energy 
density, the quintessence scenario loses its appeal. In con- 
trast, if we detected an evolution of the dark energy density, 
quintessence would likely become a very serious candidate. 



2.5.3. Alterations to gravity 

Since dark energy is required when assuming that GR holds 
on large scales, one may question the theory itself. In the GR 
framework, the evolution of structure formation is sensitive 
to the expansion history and hence constitutes an alternative 
probe of the Universe content. Precision measurements of the 
growth of structure in the past constitute a key test of the dark 
energy paradigm in the GR framework and a strong justifica- 
tion of dark energy space missions (see e.g. [46, 47]). Corre- 
lations of the cosmic shear (e.g. [48]) are a dark energy probe 
sensitive to both structure formation and expansion history. 

Before altering GR, one should note that it successfully 
passed demanding tests in the solar system (e.g. [49]) and 
that alterations to GR should only affect cosmological scales 3 . 
Embedding our 4D space-time into a 5D Minkowski space 
provides some room for alterations to standard GR ([53], and 
[54] for a simplified introduction). Although one could origi- 
nally fit the data without dark energy in this framework [54], 
improved data sets now rule out some specific 5D models 
[55, 56]. 

The standard 4D GR action integrates % (the scalar cur- 
vature, not the scale factor) and alterations to this action are 
called "/(??)" theories and have been considered as alterna- 
tives to dark energy. Quoting [57], "All f(K) modified gravity 
theories are conformally identical to models of quintessence 
in which matter is coupled to dark energy with a strong cou- 
pling", while in the simplest effective dark energy models 
(such as A), matter and dark energy are coupled only grav- 
itationally. One should hence expect a significantly differ- 
ent phenomenology: comparing predictions to data, several 



2 This 10 120 ratio of energy densities comes from summing the modes of 
quantum oscillators up to the Planck scale. This is arguable, but setting the 
cutoff to the electroweak scale "only" reduces the discrepancy to ~ 10 60 . 
In order to reconcile the zero point energies with cosmological densities, 
the cutoff should lie in the meV regime [37]. 



3 The "Pioneer anomaly" [50, 51] is the only evidence a GR "failure" on 
small scales. Improved models of radiative transfers on the spacecraft have 
recently been presented which may account for the whole effect [52]. 
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works [57-59] exclude a wide class of fCR) models and dis- 
favour those as an alternative to dark energy. 



2.5.4. Back-reaction 

The standard derivation of Friedman equations assumes a 
homogeneous Universe which is obviously a simplification. 
In [60], it is shown that in an inhomogeneous Universe, a spa- 
tially averaged Friedman equation may convey a term causing 
acceleration, only due to matter inhomogenei ties. This "back- 
reaction" mechanism relies on the non-linear structure of GR, 
and at variance with all other dark energy incarnations does 
not involve new physics. In this scheme, dark energy is a "mi- 
rage" from structure formation and the "coincidence problem" 
vanishes. 

Since one should question whether alterations of the av- 
erage Friedman equation indeed alter the average distance- 
redshift relation in the same way as for homogeneous mod- 
els, explicit toy models have been built. By "punching" low- 
density spherical voids in a homogeneous matter background, 
with a continuous density, one obtains a Swiss-cheese Uni- 
verse, for which analytic solutions have been derived (see e.g. 
[61] and references therein). In a periodic lattice of spheri- 
cal voids, distances differ significantly from a homogeneous 
dark energy free Universe, but do not mimic the concordance 
model distance-redshift relation (figure 5 of [61]), even with 
voids as large as 350 Mpc. Randomizing void locations makes 
the disagreement even worse [62]. 

Another stream of back-reaction attempts develops second- 
order perturbation theory in order to evaluate the contribu- 
tion of perturbations to source terms of Einstein's equations 
(see e.g. [63, 64]). Although one finds contributions that 
would contribute as some effective dark energy in the Fried- 
man equations, their physical meaning and gauge invariance 
have been questioned [65]. It has also been argued that correc- 
tions to large-scale distances from density contrasts are very 
small (e.g. [66]). 

Because back-reaction is technically very involved, 
whether it constitutes an alternative to dark energy is debated: 
compare, e.g., [67] and [65]. One should, however, remark 
that (perhaps for excellent reasons, [67]) the back-reaction 
mechanism does not deliver yet the ACDM phenomenology 
without dark energy. 

3. ASTRONOMICAL OBSERVATION TECHNIQUES 

This short section, for readers unfamiliar with astronomy, 
introduces the minimal concepts and parlance of astronomical 
observation at optical wavelengths, and we refer to [68] for a 
textbook devoted to observation techniques. 

Astronomical observations at optical wavelengths consist 
mostly in imaging and spectroscopy. Imaging is amost always 
carried out through filters, in particular for measuring dis- 
tances to SNe la. There are two sets of common filters in the 
visible, the UBVRI "Standard filters" [69] and the ugriz SDSS 
filters [70], displayed in figure 1. Practical implementations 
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FIG. 1. Shape of commonly used photometric bands in the visible, 
arbitrarily normalised. U, B, V, R, I refer to standard bands from 
[69], and u, g, r, i, z to the SDSS bands from [70]. The latter integrate 
representative atmospheric absorptions experienced from the ground. 



of filters are always slightly different from the archetypes. In 
the visible range, light sensors are now essentially always sil- 
icon CCDs, which have a red cutoff at around 1.1 /im and 
deliver pixelized images that can now reach millions of pix- 
els per device. Before the 1990s, images were mostly ac- 
quired with photographic plates (in particular wide-field im- 
ages), and precision photometry used vacuum tubes, typically 
photo-multipliers. CCDs became common during the 1990s 
and wide-field imaging is now carried out using mosaics of 
CCDs, with tens to thousands of individual devices. Typi- 
cal imaging observations consist in several minutes of inte- 
gration (where the integration time is defined by a mechanical 
shutter) followed by about 1 min of read out. Typical read- 
out cadences are O(100 000) pixels per second, in order to 
limit the contribution of electronic noise to ~5 electrons per 
pixel. The dominant noise affecting ground-based CCD im- 
ages is most usually the shot noise of the sky glow, except 
for the brightest sources where the shot noise of the object 
itself contributes. Image quality (IQ), which directly deter- 
mines the photometric accuracy for point-like objects such as 
stars (and supernovae), is usually described through the size of 
these point-like sources, and the best ground-based sites rou- 
tinely reach IQ below 1" FWHM in the visible. This figure is 
almost one order of magnitude better in space. Astronomers 
usually report flux measurements as magnitudes defined as 
m = -2.5/ogio(flux) + C, where C is defined to attribute a 
certain magnitude to some reference object, either real (such 
as Vega) or some spectrum such as AB magnitudes (see [71]). 
How these magnitude systems are related to physical fluxes is 
described in §9. "Colour" refers to magnitude combinations 
insensitive to overall brightness, such as B—V. For quasi-static 
stars, colours are a rough proxy for temperature. 

For spectroscopy, one relies on a dispersing element, com- 
monly a grating. To select the light from the object of interest, 
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the most common technique consists in positioning a narrow 
slit in the image plane and then dispersing the light in the di- 
rection perpendicular to the slit. The slit can be replaced by 
fibre optics [72] or systems of lenslets [73]. Multi-object spec- 
troscopy (MOS) consists in simultaneously acquiring spectra, 
either using fiber optics positioned on the targets in the tele- 
scope focal plane, or machined multi-slit masks. Currently, 
high-multiplex MOS instruments typically acquire 0(1000) 
spectra at a time. In space, one can perform slit-less spec- 
troscopy by placing a disperser in front of an imager, and 
still separate the object light from the sky light, thanks to the 
low sky brightness. The spectral resolution 5A is commonly 
expressed through R = A/ 6 A. As imaging, spectroscopy in 
the visible also uses CCDs. Synthetic photometry from slit 
spectra is generally inaccurate, because the fraction of light 
missing the slit is both unknown and wavelength dependent. 
Some spectrograph designs explicitly tackle this limitation 
(e.g. [73, 74]). 

Observations in the near infrared (NIR) now rely on low 
band-gap pixelized semiconductor devices (e.g. [75]), cou- 
pled to integrated readout electronics. Observing around and 
above 1 //m from the ground is difficult, because the atmo- 
sphere glow rises very rapidly with wavelength, and also fea- 
tures numerous absorption lines. Both atmospheric glow and 
absorption are variable on timescales of minutes to hours. 

Observations from space benefit from a much improved IQ. 
Compared to ground, the sky background is reduced, with a 
dramatic improvement in the NIR [76]. On the other hand, 
ionizing radiation hitting the sensors affects the acquired data 
via local charge deposition, and also degrades the perfor- 
mance with integrated dose. 



4. PIONEERING WORKS ON SUPERNOVAE AND THE 
PATH TO CURRENT CLASSIFICATION 

Already in the late 20's, Baade and Zwicky were search- 
ing transients by visually comparing photographic plate pairs 
exposed for about 30 min, about one month apart, on the Palo- 
mar 18 inch telescope, in order to identify transient events. In 
1934, they proposed the name "supernovae" [77] to denote a 
new class of extragalactic "novae" typically 10 4 times brighter 
than regular novae. They declared these two event classes as 
totally different phenomena and proposed a supernova rate of 
about one event per galaxy and per several centuries. They 
identified that nuclear explosions are necessary to generate the 
enormous amount of energy they observed. They were able to 
acquire spectra of their events, which exhibit the broad fea- 
tures characteristic of fast moving ejecta, and which do not 
show evidence of hydrogen lines. In 1941, a second class of 
supernovae was proposed [78]: these "type II supernovae" ex- 
hibit hydrogen lines and are typically fainter than the original 
"type I" supernovae. The classification of type I supernovae 
was refined in 1985 [79], when a subclass of peculiar type I 
were isolated: those are typically fainter than regular type I 
and miss the Si II broad absorption at ~610 nm. Type la su- 
pernovae were then defined as events that display Silicon and 
no hydrogen whereas type lb and and Ic do not display silicon 



nor hydrogen. Type lb display a strong helium line, and type 
Ic do not display it. 

The current taxanomy of supernovae was sketched in 1993 
[80], and its current status and history are detailed in [81]. 
From a physical point of view, type la supernovae are ther- 
monuclear explosions, while all the other types are core- 
collapse explosions. We refer to [82, 83] for introductory and 
more advanced material on theory of supernova explosions. 

Core-collapse progenitors are massive stars (typically 
above 8 M [84]) which run "out of fuel" when the chain of 
successive fusions stops on Fe nuclei: further fusions are en- 
dothermic. The gravitational pressure can then no longer be 
balanced by heat radiation, and the core collapses (e.g. [85] 
and references therein). Iron nuclei are photodisintegrated and 
protons are turned to neutrons by electron capture and neu- 
trino emission, producing a neutron star, on which the outer 
shells fall and bounce. The spectral characteristics depend 
mainly on the composition of the outer shells, where lb and Ic 
have been stripped of H and/or He, and type II have retained 
their hydrogen. Since the progenitors are massive, they expe- 
rience a rapid evolution (lifetimes of 10 6 - 10 7 years), so that 
the core-collapse SN rate traces the star formation rate. Core- 
collapse supernovae arise from a variety of progenitors and 
exhibit a broad distribution of luminosity, even within each 
subclass [86, 87]. As already noted, these events are on av- 
erage fainter than SNe la [87], and the population of the faint 
tail is still poorly known. 



5. OBSERVATIONAL CHARACTERISTICS AND 
DIVERSITY OF SNE IA 

5.1. Light curves and colours 

SNe la emit most of their energy in the visible, and the 
remaining part is essentially confined to near UV and near IR. 
During the bright phases, the spectrum peaks around 400 nm 
(see Fig 3). Very rapidly, observers noted that SNe la (which 
were still called SNe I) exhibit reproducible light curves (e.g. 
[88]). 

In the B-band (blue light, see figure 1), the light curve 
rises in about 19 (restframe) days on average (e.g. [89]). 
About 35 days past maximum the light curve follows a quasi- 
exponential decay (see figure 2). If one sums the emitted en- 
ergy, thus building the bolometric light curve, the observed 
decay rate beyond ~ 50 days after maximum is related to the 
56 Co half-life: the B and y decay chain 56 Ni ~^ 56 Co -^ 5b Fe 
powers the light curve (e.g. [90-93], and §6). 

In the /?-band (~650 nm), the light curves exhibit a shoulder 
20 to 30 days after maximum, which transforms into a second 
maximum towards redder bands (up to ~2fim), taking place 
20-30 days after the first one (see Fig 2); the feature is more 
prominent for bright events, faint events might miss it. 

The B - V colour (blue versus green, where V stands for 
"visible", see also figure 1) reddens during the first 50 days, 
indicating that the emitting material is cooling. At about 40 
days after maximum (i.e. ~ 60 days after explosion) the B-V 
colour changes slope and then becomes bluer and bluer (e.g. 
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figure 4 of ([94], figure 1 in [95]). Other colours have different 
behaviour. 



5.2. Spectra 

SNe la spectra captured early or around maximum are char- 
acterized by broad absorption lines of "intermediate-mass" el- 
ements, mainly Si and Ca, shifted to the blue with respect 
to the SN host galaxy redshift, because the absorption takes 
place in layers ejected at high velocity towards the observer. 
The velocity of the absorbing region is measured from the 
blue-shift of absorption features in the spectrum and decreases 
with time (e.g. [96, 97], and figure 3). The light-emitting re- 
gion (called photosphere) is characterized by a material den- 
sity below which photons can escape. As the supernova ex- 
pands, its density decreases and the light-emitting region pro- 
ceeds towards the core of the star that just exploded (its nature 
will be discussed at §6). The supernova eventually enters a 
nebular phase where emission features dominate. 



5.3. Variability 

SNe la exhibit a variability of the width of their light 
curve, even after accounting for time dilation in order to 
compare restframe widths. This variability was initially no- 
ticed through a correlation between the decay slope of the B 
(~450 nm) light curve and its peak brightness ([100, 101]). 
This correlation was refined in 1993 from a small sample 
of well-measured events [102]. The "brighter-slower", or 
Phillips relation quantifies that slower decliners are also in- 
trinsically brighter, which is of paramount importance when 
measuring distances: all supernova cosmology analyses im- 
plement it in some way. Still in the B-band, it has been pro- 
posed to describe this variability of light curves by stretching 
the time axis of a single light curve template [103, 104], and 
the rise and fall timescales seem to vary together for [104] and 
[89], while [105] finds them essentially independent. 

The distribution of decline rate (or stretch factor) has an rms 
variation of ~10% [106], exhibits a long tail of hard-to-catch 
rapidly declining and faint events [107], as well as some bright 
and slow events. Cosmological analyses have to exclude 
these extreme events, because the empirical light curve shape 
parametrizations behave poorly for those [106, 108, 109], and 
also because the brighter-slower relation changes [107, 110]. 
Fast decliners are on average redder than slow decliners [1 10— 
112]. 

We will discuss parametrizations of this light curve shape 
variability when we come to light curve models in §8.3 but all 
turn out to be essentially equivalent, and all extract the relation 
of their parameter to luminosity from data. 

SNe la also exhibit variability of their colours (measured, 
e.g., at maximum) even at a fixed decline rate. For exam- 
ple, the distribution of the restframe B-V colour has an rms of 
about 0.1 mag. (see e.g. [106]) and a tail of highly reddened 
events (e.g. figure 5 in [113]). Different colours of the same 
event seem to be well correlated: B-V and U - B are related 



by a linear relation (see e.g. figure 7 in [106]), as expected if 
colour variability is attributed to varying amount of absorption 
along the line of sight. More precisely, colour relations of SNe 
(such as U-B versus B-V) slightly depend on the decline rate, 
which is easily accounted for (e.g. [109]). As one considers 
bands bluer than B, the relation between different colours be- 
comes poorer and poorer (e.g. figure 7 in [114], [115, 116]), 
a tendency confirmed by space-based UV photometry down 
to 180 nm [117]. As for extinction, one observes that bluer 
supernovae are brighter, and that different colours do not vary 
independently: one colour seems enough to parametrize this 
colour variability. The source of colour variability unrelated 
to the decline rate is still debated: it could be intrinsic to su- 
pernovae (e.g. [118]), it could be due to extinction by dust 
in the host galaxy, or more likely a mixture of both. There 
is evidence of narrow Na absorption lines, indicative of dust, 
particularly for highly reddened events (e.g. [113] and refer- 
ences therein), which are usually excluded from cosmologi- 
cal analyses. One might, however, note that the quantitative 
agreement of SNe colour variations with known dust proper- 
ties is poor (see e.g. [108, 109, 112, 119]), with one notable 
exception [120]. We will discuss relevant parametrizations in 
§8.2 & 8.3 . 

The brighter-slower relation becomes shallower when go- 
ing to redder wavelengths [121], and eventually vanishes at 
NIR wavelengths [122, 123]. The brighter-bluer relation 
also becomes shallower, and ignoring both of them deliv- 
ers accurate distances in NIR bands in particular in the ri- 
band ([1.5,1.8]/mi) [123] (but see also [124] and references 
therein). These findings are so far limited to nearby events, 
because these NIR restframe bands become increasingly diffi- 
cult to observe as redshift increases. 

SNe la variability is often described as a one-parameter 
family, with extinction by dust in the host galaxy causing 
some extra colour variation. In practice, and whatever the 
source of colour variations, cosmological works mostly con- 
sider these events as a two-parameter family, where one pa- 
rameter indexes light curve shape and the other some colour, 
or colour excess. 



5.4. Correlations with host galaxy properties 

The physics in galaxies is broadly indexed by the mean age 
of their stellar population. Early-type galaxies mainly con- 
tain evolved stars and form little or no new stars, while late- 
type galaxies have varying amounts of star formation going 
on. One might broadly rank galaxies from big, red, passive 
and elliptical to small, blue, star-forming and spiral or irreg- 
ular. SNe la are observed in all galaxy types, while core- 
collapse supernovae are observed only in star-forming hosts 
(see, however, [125]) because their massive progenitors have 
a very short lifetime. 

The SNe la rate can be split across the galaxy types (see 
[126, 127] and references therein), and one finds that star- 
forming galaxies produce more than 10 times more SNe la 
(per unit stellar mass) than passive galaxies. Fast decliners 
happen preferably in old stellar environments, or early-type 
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FIG. 3. Spectral sequence of SN2003du ([99]). Spectra are labelled 
by their phase w.r.t B maximum. Absorption features are attributed 
to Si and Ca elements. Note that the absorption features of the first 
spectra shift to the red with time. 



hosts (figure 4. of [9], e.g. [127-129]), and the cause remains 
unclear. This probably explains the early finding that SNe la 
are brighter in late-type galaxies [130], and the observation 
that the average decline rate decreases with redshift [131]. 
Colours have been compared across environments and found 
to be compatible [128, 132], which is somehow surprising be- 
cause if extinction by dust contributes significantly to colour 
variations, one expects less extinction in passive galaxies be- 
cause they have less dust than active ones (e.g. [133, 134]). 
Even more surprisingly, the largest homogeneous SNe la sam- 
ple to date recently indicated [135] that supernovae in passive 
hosts are slightly but significantly redder than those in star- 
forming hosts. A similar trend was reported earlier in §4.2.2 
of [136]. Some other aspects related to host galaxies are dis- 
cussed in §10.1.3. 



6. THE PARADIGM OF SNE IA EXPLOSIONS 

6.1. Progenitors and initial state at ignition 

The current model for type la supernovae explosions as- 
sumes that a carbon-oxygen (C-O) white dwarf (WD) accretes 
mass from a companion until it reaches the Chandrasekhar 
mass {Mchcm = 1.38M©), at which the pressure from the 
degenerate electron gas no longer balances the gravitational 
pressure [137]. The temperature can then increase enough to 
trigger carbon fusions [138], and the nuclear flame can "burn" 
the whole star material within a few seconds. A large con- 
sensus has built up around this paradigm [139], because it is 
supported by several observational facts. 

• Fusions of C and O nuclei produce Si, S and Ca, which 
are observed in the spectra (see figure 3). 



• Ni is the expected end of a fusion chain starting with 
C and O since further fusions would not release energy. 
The 56 M(T l/2 = 6.1d) -^ 56 Co(T 1/2 = lid) -> 56 Fe 
seems to power the light curve, both because of its de- 
cay time (e.g. [90-93]), and because of the evolution of 
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line strength in the late spectra [140]. 

• The reproducibility of luminosity (e.g. [121]) is a logi- 
cal consequence of the explosion occurring at a thresh- 
old depending only on fundamental constants. 

• The lack of H and He in the spectrum is just due to their 
absence in the initial state. 

• The fact that SNe la are the only type of SNe observed 
in old stellar environments such as elliptical galaxies 
suggests a low-mass (hence long-lived) evolved progen- 
itor. White dwarfs form at the end of the evolution of 
stars of less than 8 M , which eject most of their outer 
layers, leaving a ~ IM G C-O compact star, too cool to 
host nuclear reactions, and slowly radiating its energy, 
possibly for billions of years. 

This scenario is strongly supported by a quantitative study 
of the outcome of nuclear reactions [17, 141], It, however, 
faces a major issue: the path to the initial state is still ambigu- 
ous ([18, 139] for detailed discussions of the various issues). 
Since the Chandrasekhar mass has to be reached from below, 
a mechanism for the WD to gain mass should be proposed. 
Binary systems (i.e. bound systems of two stars) provide two 
different paths, the so-called single-degenerate (one WD and 
some other star) and double-degenerate (two WDs) scenarios. 
A recent survey of WDs in Sun's neighbourhood [142] reports 
that 25% of those are in binary systems, and 6% in double- 
degenerate systems. Among known low-mass WD binaries, 
more than 10 systems will merge in less than a Hubble time 
[143]. 

The single-degenerate scenario consists in a low-mass WD 
accreting material from a higher mass companion star until 
it reaches the Chandrasekhar mass (see e.g. [144]). This re- 
quires that the mass donor is a low-surface-gravity and non- 
degenerate star, and late main-sequence stars or red giants are 
the most likely candidates. H and He are the more likely ele- 
ments to be accreted. However, plausible accretion rates face 
serious issues (see p 206 of [17] and reference therein): slow 
accretion rates do not compensate for the mass loss due to no- 
vae eruptions, or helium flashes; in contrast, at high accretion 
rates, the steady burning of H and/or He should make these 
progenitors so bright that they should be easily detectable. 
The low observed X-ray flux from elliptical galaxies might 
forbid [145], or not [146], the single-degenerate scenario, at 
least in elliptical galaxies. However, the arguments relying on 
the luminosity induced by material accretion may not apply 
if the progenitor binary systems are enshrouded during their 
accretion phase. We should stress here that the binary sys- 
tem (and its environment) is complicated enough to take the 
outcome of models with caution. For example, the role of H 
and/or He flashes is also controversial [147, 148]. The single- 
degenerate scenario is more and more challenged by the donor 
star escaping observations: none was found in pre-explosion 
archival images [149, 150], nor within SN la remnants (e.g. 
[151] and references therein); the expected UV and X-rays 
emission from the explosion material colliding with it was 
not seen [152, 153]. We cannot entirely rule out the single- 
degenerate scenario, but its viability will be questioned if the 



donor star remains elusive. Note that WD are generally far too 
faint to be detected in pre-explosion archival images. 

The double-degenerate scenario assumes an initial binary 
system of two WDs [154, 155], losing orbital angular momen- 
tum by means of gravitational radiation. Stellar population 
synthesis calculations predict the required statistics for merg- 
ers [139]. While this scenario naturally explains the absence 
of hydrogen in the spectrum, there is no built-in mechanism 
that leads to a reproducible state at ignition, a simple expla- 
nation for the empirical reproducibility of SNe la. On the 
other hand, this initial state naturally provides a mechanism 
to explain extremely luminous events such as SNLS-03D3bb, 
which seems to require a super-Chandrasekhar initial state to 
produce enough 56 M to explain its apparent luminosity [156]. 
However, the single-degenerate initial state might also gener- 
ate supra-Chandrasekhar amounts of 56 Ni [157]. A threat to 
this double-degenerate initial state is that WD mergers lead to 
off-centre ignition causing gravitational collapse rather than 
the expected thermonuclear disruption (e.g. [83, 158, 159]), 
but a recent simulation actually obtains an under-luminous 
SNe la from two 0.9 M WDs [160], but requires a narrow 
mass difference of the involved stars. 



6.2. Explosion models 

Assuming the initial state can indeed be obtained, mod- 
elling the outcome of the explosion is at least a two-step pro- 
cess: one first has to go through the nuclear fusion reactions, 
and second to evaluate the radiation transfers in the expand- 
ing ejecta. Both aspects are physically and computationally 
extremely challenging [17]. 

The production of elements through successive fusions dur- 
ing the thermonuclear runaway depends on the reaction con- 
ditions, mainly density and temperature: the reaction rates 
vary extremely rapidly with temperature, and the production 
of iron-peak elements 4 favoured under high density. There are 
roughly two regimes for the propagation of the flame: defla- 
gration (subsonic) and detonation (supersonic), where defla- 
gration allows the unburnt material to expand before reacting. 
Detonation models (see e.g. [161]) very efficiently convert the 
input material into 56 M, and then under-produce intermediate- 
mass elements (Si, Ca, S,...) which constitute the prominent 
features of the spectrum around maximum light. The famous 
W7 explosion model [162] produces intermediate-mass ele- 
ments via deflagration, and is a fair match to the observed light 
curve and spectrum at maximum. So, a deflagration phase is 
required in order to reduce the density of the stellar material 
prior to nuclear reactions in a Mc/ M „ WD [17, 163], in order to 
avoid the nearly exclusive production of 56 M [83, 162, 164]. 

However, pure deflagration models face some discrepancies 
with observations [162, 163]: overproduction of 54 Fe w.r.t to 
solar abundances (sometimes referred to as over-production 



"Iron peak" refers to a local maximum, around Fe, in the abundance of 
elements as a function of atomic number. The expression refers here to 
Fe.Co and Ni. 
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of neutron-rich elements), some indication of a brighter-faster 
relation (while we observe a brighter-slower relation), and too 
narrow an interval of expansion velocities of intermediate- 
mass elements. To cure these issues, [163] proposes that af- 
ter an initial episode of deflagration, the flame transitions to 
detonation (deflagration-detonation transition (DDT)), as ob- 
served in laboratory chemical combustion experiments (e.g. 
[165]). Through the alliance of both combustion modes, 
intermediate-mass elements are produced over a large range 
of ejection velocities. DDT simulations hence improve the 
match to observed light curves and spectra [166], and provide 
some hints that the density at transition might source the ex- 
plosion strength variability [166, 167]. 



6.3. Radiative transfer 

The explosion determines the chemical composition of the 
ejecta and velocities of the various layers, mostly driven by 
the energy released in the nuclear fusions. From this interme- 
diate state, radiative transfer calculations generate synthetic 
light curves and spectra, which are mostly powered by the 
56M /3 and y decay chain (e.g. [90-93]) as already dis- 
cussed. Broadly speaking, the supernova emits less energy 
than it produces until the maximum, when both quantities are 
equal [161]. Arnett's rule ([168]), derived from analytic ra- 
diative transfer analysis assuming a grey opacity, states that 
the maximum light luminosity is proportional to 56 M mass 
and has been used to constrain the latter (e.g. [169]): it varies 
by about one order of magnitude from very fast to very slow 
decliners, and by about a factor of 2 for the range of decline 
rates commonly included in cosmological analyses. 

Radiation transfer codes simulate, in an expanding and in- 
homogeneous medium, the conversion of nuclear radiation 
into the eventually observed light; their sophistication has fol- 
lowed the increase in available computing power. The com- 
plexity arises from two causes: first, opacities depend on the 
ionization states of the atoms, which themselves depend on 
the radiation flux and temperature, which in turn depend on 
opacities. One thus has to find a consistent solution, and the 
convergence speed is a major concern (see [170] and refer- 
ences therein). Second, the opacities are driven by interactions 
with atomic lines (scattering and fluorescence), which are ex- 
tremely numerous. In particular, the iron group elements ex- 
hibit a wealth of lines in the near U V that mostly shift light to 
the red, but also produce most of the UV flux around 300 nm 
([171, 172]). Both the details of "atomic data" (i.e. the choice 
of lines involved in the calculation) and the treatment of ion- 
ization are shown to have a large influence on the predicted 
light curves and spectra (e.g. figure 7 of [173]). 

The current state of the art mostly consists in 3D solu- 
tions. Monte Carlo techniques [174] are efficiently paral- 
lelized and the microscopic physics can be incrementally re- 
fined. The various implementations (e.g. [173, 175] and refer- 
ences therein) mostly differ in the details of the micro-physics 
approximations. Hydrodynamic solutions now also accomo- 
date 3D (e.g. [176] and references therein). Going to three 
dimensions allows one to account for possible anisotropics, 



which have been proposed as a potential source of diversity 
through the random viewing angle [173, 177]. 



6.4. Models at late times 

At about 100 days after maximum light, emission lines 
start to dominate the spectrum, indicating that the luminous 
material is scintillating because of the energy deposition of 
positrons and y rays. This feature has been used to constrain 
both the amount of ionizing radiation and the nature of ma- 
terial around it, using models much simpler than in the early 
phases. At t>~200 days, the quasi-exponential slope of light 
curves is different from band to band [178, 179], possibly due 
to the increasing amount of (stable) 56 Fe which scintillates 
significantly in the B-band [178], and also because more and 
more y-rays and positrons escape [91, 180]. The positron es- 
cape fraction constitutes a handle on the magnetic field [92]. 
Simulations of late SNe la [180] and y-ray observations [181] 
indicate that "late" positrons from SNe la might constitute a 
dominant contribution to the 511 ke V emission in the Galaxy. 
From the difference of the bolometric late light curve slope 
with the 56 Co decay rate, one can model the escape rate of y 
rays and obtain constraints on the ejected mass (e.g. [91, 93]) 
and on the amount of synthesized 56 M ([169] and references 
therein). 



6.5. Explosion models and cosmological use of supernovae 

The value of Ho has been inferred from a Hubble diagram 
of nearby supernovae and numerical simulations to derive the 
intrinsic luminosity of events [166, 182, 183], and values fully 
compatible with conventional determinations (mostly relying 
on the cosmic distance ladder [19]). This approach is proba- 
bly the sole quantitative contribution to date of SN explosion 
models for determining cosmological parameters. Because 
of the complexity of the explosion and the consequent light 
production, models nowadays only broadly reproduce light 
curves and spectra and cannot be used as templates for fitting 
data. If the brighter-slower relation is qualitatively understood 
and physically attributed to synthetized 56 M mass variations 
(e.g. [184] and references therein), there is no consensus on 
what causes the diversity of the latter (e.g. [185]). Models 
hence do not yet provide strong guidance to implement the 
brighter-slower relation. The brighter-bluer relation is a worse 
case because there is no consensus on whether variability of 
supernovae (independently of the decline rate) contributes sig- 
nificantly to it (see §8.2). 

Cosmologists, however, hope to get some insight from fur- 
ther improving the physical understanding of the initial state, 
its environment and the explosion mechanism(s). Ideally, we 
would improve light curve modelling, reduce the distance 
scatter, and gain insights into redshift-dependent systematic 
biases of distances. Such goals might look excessively ambi- 
tious, but we should remark that qualitative indications might 
be sufficient. For example, if models point out some measur- 
able indicator of the WD metallicity, the relation of metallicity 
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to intrinsic luminosity can be inferred from data. 



7. OBSERVING AND REDUCTION TECHNIQUES FOR SN 
SURVEYS 

Up to now, most SNe la events that were used to measure 
distances went through three different observing phases: find- 
ing the event, identifying it from a spectrum (spectroscopic 
follow-up), and measuring its light curve (photometric follow- 
up), possibly in several bands. We will now discuss how the 
required observations are conducted and handled, with some 
historical perspective. 



7.1. Finding supernovae 

Essentially all supernovae are nowadays found using an ap- 
parently simple technique: image subtraction. The concept 
consists in subtracting a search image from a previously ac- 
quired image often called reference. The subtraction is done 
pixel to pixel and point-like positive excesses on the resulting 
image constitute supernova candidates. If the time gap be- 
tween the two epochs is chosen to be around 2 weeks (in the 
restframe), found supernovae are likely to be still rising and 
a prompt follow-up campaign allows one to capture the max- 
imum light (which allows one to safely estimate peak bright- 
ness, and eases spectroscopic identification). Prior to the sub- 
traction itself, both images have to be brought on the same 
pixel grid, on the same flux scale and their point spread func- 
tion (PSF, in practice the shape of stars, which varies with 
atmospheric conditions) should be homogenized. The pro- 
cess aims at making non-variable objects disappear in the sub- 
traction process. Since the concept of image subtraction was 
initially proposed [186] and demonstrated on the sky [187], 
the major improvement brought to the technique consists in 
fitting the convolution kernel for PSF homogenization [188], 
by explicitly minimizing the subtraction image. Subtractions 
are vulnerable to small defects of geometrical alignment or 
PSF matching, or more generally, significantly different ob- 
serving conditions, and reaching the shot noise limit requires 
some care. Image subtractions are then essentially always per- 
formed using image pairs from the same instrument. One has 
to promptly find candidates in order to rapidly obtain spec- 
troscopic confirmations and start the photometric follow-up 
as soon as possible. The computing load involved in find- 
ing supernovae has not decreased significantly in the two two 
decades since the area of silicon sensors (nowadays often as- 
sembled into mosaics) used to image the sky, and the comput- 
ing power of silicon processors used to find supernovae have 
followed similar slopes. An example of two images tiles and 
their subtraction is displayed in figure4. 

Searching supernovae requires an imager with as large a 
field of view as possible. The first sizable SNe la sample (at 
z<0.1) with CCD-measured light curves [9] was found us- 
ing photographic plates in the early 1990s. Discovering su- 
pernovae in a somehow systematical way using CCDs started 
around 1995: for example, the IAU telegram 6270 announced 
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FIG. 4. Small portions of images of the same field. The leftmost 
was acquired on 8 September 1999 and gathers 3.6 h of integration 
on CFHT, using the CFH12k imager in I band. The central tile was 
taken 26 days later with the same instrument. The rightmost tile dis- 
plays the image subtraction, where a positive excess is visible at the 
centre, at a brightness I Vega ^ 24.7 and about 6 cr. This excess was 
spectroscopically identified as a type la event (SN1999fd) at z=0.87. 



11 supernovae at 0.35 < z < 0.65, found using a 16'xl6' 
single-CCD imager on the CTIO-4m, and spectroscopically 
confirmed at the Keck (10 m). We note here that a sizable 
fraction of the distant events reported in the accelerated expan- 
sion papers [10, 11] (§11.2) was discovered using this CTIO- 
4m/Keck combination. 

The SN finding techniques and instruments then improved 
rapidly: in 1999, the IAU telegram 7312 reports 20 spectro- 
scopically confirmed high-redshift supernovae (among which 
3 events at z > 1), discovered using a 30'x40' imager at the 
Canada-France-Hawaii Telescope (CFHT, 3.6 m), and con- 
firmed at Keck. Unfortunately, the photometric follow-up of 
this wonderful harvest could not be fully conducted, illustrat- 
ing that actually measuring light curves of distant supernovae 
was becoming the heaviest part of the observing programme. 



7.2. Rolling searches 

In this 1999 discovery telegram (IAU 7312), one can eas- 
ily check from the event coordinates that several supernovae 
were found on average in the same telescope pointing. This 
suggests the "rolling search" approach, which consists in re- 
peatedly imaging the same sky patches, and using the im- 
age sequence both to find supernovae and to measure their 
light curve. The spectroscopic identification phase remains 
unchanged with respect to the "traditional" two-epoch search. 

The technique becomes increasingly efficient as more su- 
pernovae are measured in the same telescope exposure on 
average. In order to benefit from this multiplex advantage, 
three ground-based surveys implemented this technique (see 
§11.5): ESSENCE on the CTIO-4m with the 36'x36' Mo- 
saic camera (2002-2008), SNLS on CFHT (3.6 m) using the 1 
deg 2 Megacam imager (2003-2008) and the SDSS SN search 
(2005-2007) on the SDSS 2.5 m with its 1.5 deg 2 fast scan- 
ning imager. These three projects were successful in find- 
ing and measuring events, essentially at the expected rate. 
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Rolling searches using the Hubble Space Telescope (HST) 
imagers have also been implemented (§11.4), however, with 
a lower multiplex advantage than the ground-based wide-field 
imagers. 

The rolling search technique not only efficiently addresses 
the observing time shortage that rules astronomic observation, 
but also improves the data quality over the classical two-epoch 
search scheme followed by photometry. First, SN light curves 
are sampled independently of their phase, since observations 
are not conducted for a specific event; in particular, SN light 
curves contain very early points, taken even before the event 
was actually detected. Second, the search is indeed easier than 
in the early times because photometry for distance estimation 
requires deeper images than for mere search. One also usu- 
ally has more than a single image and epoch to turn a light 
excess into a candidate. Finally, these repeated images can be 
stacked to reach considerable depths, especially if the survey 
lasts several years. The wide-field stacks usually deliver the 
photometric characterization of the supernova galaxy hosts 
(e.g. [127, 128]). They can also be used for science topics 
unrelated to supernovae, which can even be the main drivers 
for the survey, as happened, in particular, for some heavy SN 
programmes on the HST (§1 1.4). 



7.3. Spectroscopy of candidates 

Candidates from supernovae searches have to be confirmed, 
most of the time using spectroscopy. From reported discover- 
ies, one can draw the crude rule that the telescope used for 
spectroscopy is twice as large in diameter as the one used 
for the search via imaging. SNe la supernovae at z < 1 are 
now found using 4 m class telescopes, and the spectroscopic 
follow-up is conducted on 8 m class facilities, where a target 
at 0.8 < z < 1 requires integrating for 1-2 h (see e.g. table 1 
in [189]). For nearby supernovae at z < 0.1, searches are typi- 
cally conducted with 1 m telescopes, and spectroscopy is done 
in 10-20 mn exposures on 2 m telescopes. In order to use the 
observing time as efficiently as possible, spectroscopic obser- 
vations aiming at classification are usually conducted around 
maximum light, if possible. For distances to SNe la, besides 
positive identification of the event against other SN types, 
spectroscopy is expected to provide the redshift value, in order 
to place the event on the Hubble diagram. 

The needed spectral resolution is modest: supernovae have 
no narrow lines, and a resolution (A/6A) of about 100 is 
enough to sample properly the peaks and troughs of SNe la 
spectra. One, however, usually uses spectral resolutions 
around 1000, in order to reduce the impact of telluric absorp- 
tions and sky line emissions, and in order to properly measure 
narrow lines of the host galaxy spectrum, which constitute the 
best redshift indicator. 

At variance with imaging, there is currently no practical 
way to efficiently multiplex spectroscopic observations of su- 
pernovae because the field of view of most spectroscopic in- 
struments is too small. The events are then observed one at 
a time. When using a multi-object spectrograph, one can still 
observe other targets (typically galaxies) in parallel with the 



SN candidate. 



7.3.1. Host galaxy subtraction 

Supernovae are point sources while their host galaxy is not. 
The size of the supernova does not depend on redshift while 
galaxies get smaller as redshift increases. Hence, the fraction 
of host galaxy light mixed with the supernova increases with 
redshift. For photometry, this is addressed by subtraction (or 
some equivalent technique, see §7.4), but for slit spectroscopy, 
subtracting a supernova-free observed galaxy spectrum is not 
practical. Moreover, most surveys do not allocate observing 
time to collect such spectra. 

Increasingly sophisticated approaches to address the host 
galactic contamination of supernova spectra have been devel- 
oped, mostly for high-redshift surveys. One can simply ignore 
the problem, thus selecting against SNe on bright hosts. One 
can attempt to synthesize the observed spectrum from tem- 
plate libraries of both supernova and galaxies [189]. In [190], 
both continuum and narrow lines are filtered out (in wave- 
length space) from the acquired spectrum, thus efficiently re- 
moving most of the contaminating host galaxy features, how- 
ever at the cost of altering the SN spectrum. In long-slit spec- 
troscopy, one can attempt to separate the point-like compo- 
nent (the SN) from a more extended background (the host 
galaxy), either with some spatial filtering [191], or by rely- 
ing on the distribution of galactic light along the slit estimated 
from multi-band imaging [192]. In [115], galaxy models are 
fitted to supernova-free multi-band imaging measurements of 
the host galaxy for subsequent subtraction from the observed 
spectrum. 



7.3.2. Type and redshift 

At high S/N, the identification of a supernova from its spec- 
trum does not usually require any sophisticated method; prob- 
lematic cases occur for events departing significantly from the 
archetypes. At higher redshift (and smaller S/N), the SN iden- 
tification is carried out by comparing its spectrum to high-S/N 
archetypal templates, for various phases and SN types, where 
the quality of the match is often quantified via least squares 
(e.g. [189, 190]), possibly matching both light curve and 
spectrum [116]. When present, narrow emission or absorp- 
tion lines of the host galaxy provide the redshift or at least 
redshift hypotheses, but one may also add the redshift to the 
searched parameter space. Redshifts are typically measured 
to a precision of 0.001 from galactic features and to 0.01 from 
SNe features (e.g. [189], figure 1 in [193], figure 7 of [194]). 
The quality of the measured spectrum obviously plays a key 
role in the identification confidence level. Regarding SNe la, a 
positive identification consists in the observation of Si and/or 
Ca spectral features (see e.g. [189] for a practical definition of 
identification confidence levels), figure 5 illustrates the sepa- 
ration of galaxy and SN spectra and the obtention of redshift 
and type of a high-redshift SN candidate. 
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FIG. 5. On the left, measured spectrum (in black) of the SNLS- 
04D4ib supernova candidate on a relatively bright host, acquired with 
FORS1 on the VLT, at maximum light. The redshift z = 0.704 is 
unambiguously obtained from the 6800 A(observer frame) jump due 
to the galaxy spectrum, confirmed by two H narrow absorption lines 
at 7000 A(H-(5) and 7400 A(H-y). A simultaneous fit to the measured 
spectrum and the light curves of the event is carried out using the 
SALT2 model (§8.3.2), where the spectrum is modelled as the sum 
(red line) of a galaxy (blue line) and SN la templates. On the right, 
the SN spectrum obtained by subtracting the best fit galaxy template 
(SO, blue curve on the left) from the measured spectrum, compared 
to the fitted SN la model spectrum (red line). This spectrum displays 
a broad absorption (characteristic of SNe) around 3800A resframe, 
attributed to Call, and was identified as a SN la because the general 
shape of the spectrum leaves no ambiguity at this phase (known from 
photometry). Figure reproduced with permission from [192]. 



7.4. Measuring light curves and photometry techniques 

The collection of photometry points along the SN light 
curve can be carried out with a small field instrument provided 
there are enough stars in the images to serve as photometric 
and geometric anchors. As indicated above, surveys target- 
ing high redshifts now use the repeated imaging approach and 
hence do not need dedicated and pointed imaging observa- 
tions. 



Supernova photometry consists in measuring the flux of 
the supernova in the survey images. Since one cannot spa- 
tially separate the light of the supernova from that of its host 
galaxy, images of the field before the supernova explosion or 
at least one year after are usually required for this task. Ex- 
plicit or implicit subtraction of these supernova-free images 
from supernova-on images is generally necessary. Rolling 
searches obtain these supernova-free images within their core 
programme if they monitor the same field for at least two 
years. 

Astronomical photometry involves two measurements: one 
has to measure the source of interest and some standard stars 
with the same instrument. So, astronomical photometry in- 
deed delivers flux ratios, and we will discuss the calibration of 
the standards themselves in §9.1. For ground-based observa- 
tions, the transparency of the atmosphere has to remain suffi- 
ciently constant between both sets of observations for the cal- 
ibration process to deliver anything useful, and astronomers 
label a night as "photometric" if the ensemble of standard 
stars observed throughout the night indicates comparable at- 
mospheric extinctions. In order to salvage non-photometric 
(or uncalibrated) observations of supernovae, one can rely on 
the assumption that stars in the same images as the super- 
nova are on average non-variable. The task of SN photometry 
hence consists in delivering the flux ratio of the supernova to 
neighbouring stars, and the calibration process will have to 
deliver, once and for all, the flux ratios of these neighbour- 
ing stars to some standards. Space-based photometry can rely 
on the usually excellent long-term stability of the instruments 
and may follow a simpler calibration path (see e.g. §2.1 in 
[195]). 

For faint sources, photon shot noise constitutes, in prac- 
tice, the main source of measurement uncertainty, and one will 
seek a statistically optimal flux estimator, in order to maxi- 
mize the number of measured events within a fixed observing 
time allocation. In fact, observing times are usually tailored to 
deliver some predefined signal-to-noise ratio, assuming an op- 
timal flux estimator. Since both the supernova and surround- 
ing stars are point sources, they have the same spatial shape, 
namely the Point Spread Function (PSF), which can be mod- 
elled on each image using bright stars. One can model the 
SN observations as a time-variable point source at a fixed po- 
sition on the sky, possibly on top of a non-varying galaxy. 
Fitting such a model to the image series delivers statistically 
optimal estimators of the flux at all epochs, together with their 
covariance matrix. Practically validated implementations of a 
simultaneous fit to an image series are described in at least 
two instances [196, 197], where the latter avoids image re- 
sampling. Other approaches to SN photometry are bound to 
be statistically less or as efficient, as stated by the Cramer-Rao 
inequality. 

Statistical efficiency is usually not a serious concern for 
nearby supernovae and the subtraction of supernova-free im- 
ages from follow-up images with subsequent photometry 
proves to be sufficient for statistical uncertainties to be sub- 
dominant (see e.g. [98]). 

Measuring point sources using their shape is referred to as 
PSF photometry, and is obviously adversely affected by im- 
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perfections of the PSF model. However, flux ratios are not 
affected by imperfections of the PSF model, as long as the lat- 
ter enters in the same way into both terms of the ratio. This is 
not the case for standard PSF photometry, optimized to fight 
shot noise (e.g. [198]), which makes the ratio of faint to bright 
stars depend on the fidelity of the PSF modelling. In contrast 
[106, 196], one might use for all brightnesses the flux esti- 
mator optimal for faint stars. Faint sources are still optimally 
measured, flux ratios become independent of the PSF fidelity, 
and fluxes are linear functions of the data. One might indeed 
regard the latter as a requirement when measuring fluxes for 
luminosity distances. 

One important concern for PSF flux measurement of faint 
sources is that fitting simultaneously flux and position in- 
volves a flux bias at low signal to noise (see e.g. appendix 
B in [106]), which might have a non-negligible impact when 
supernovae become faint. The obvious remedy consists in fit- 
ting a single common position to all images and even to all 
bands. In order to limit systematic biases of flux ratios, one 
should then enforce the same constraint for measuring fluxes 
of the surrounding stars. For surveys spanning several years, 
one may have to account for proper motions of stars [197]. 

One might have noted that SN measurements involve ex- 
plicit or implicit subtractions, while star measurements do not. 
This difference might systematically bias flux ratios. For this 
reason and for mere method validation, supernova photometry 
is commonly verified by inserting artificial supernovae (either 
synthetic, or copies of real stars) into real images, in order to 
bound the bias of recovered flux ratios below dominant sys- 
tematic uncertainties [196, 197, 199]. Rolling searches allow 
one to measure pre-explosion SN fluxes, and checking that 
those average to zero indicates that galaxy fluxes are properly 
subtracted [106, 197]. 

Although astronomers usually report magnitudes, it is 
mandatory that light curves points are reported in flux space, 
because some low-S/N measurements unavoidably happen to 
be negative, and retaining those is mandatory for flux esti- 
mates to remain unbiased. For the same reason, light curve 
fitters should operate in flux space [109, 114, 200], but might 
report their results in magnitudes (as commonly done). 



8. LIGHT CURVE FITTERS, DISTANCE ESTIMATION 

8.1. K-corrections 

In order to compare supernovae, their fluxes have to be 
expressed in the same way, and in particular using the same 
restframe bands. Different supernovae are generally not mea- 
sured using the same photometric setup: first, since the quality 
of cosmological constraints improves with redshift lever arm, 
nearby and distant supernovae are measured with different 
telescopes and hence different spectral bands. Second, even if 
all supernovae were measured using the same telescope, the 
filters used for the measurement map differently in the su- 
pernova restframe when the redshift varies. We hence need 
to somehow transform the photometric measurements carried 
out with an arbitrary setup into what some well-defined instru- 



ment, not necessarily real, would have measured. 

In the astronomical literature, transforming photometric 
data obtained with some filter into what would have been mea- 
sured with some other filter is called K-correction. It is usu- 
ally presented as a magnitude correction and was originally 
designed to place galaxies on a Hubble diagram [201]; these 
galaxies were measured in similar filters but at different (and 
moderate) redshifts. The principle is simple (e.g. equation 
(B7) in [201]): correct the flux measurement f(T 2 ) by a mul- 
tiplicative factor computed using an assumed spectrum and 
the involved band passes: 



f(Ti) = f(T 2 ) 



J S(A)Ti(A)dA 
J S(A*(l+z))T 2 (A)dA 



(6) 



where T 2 is the transmission function of the filter that was 
used for the measurement, T\ is the transmission function of 
the filter that would have been used ideally, and S is the as- 
sumed restframe spectrum of the object (at redshift z). Orig- 
inally, the concept was developed in a framework where T\ 
and T 2 denote the same band. Obviously, the transformation 
depends more and more on the knowledge of S as the overlap 
between T\(A) and T 2 (A/(l + zj) decreases. Note that the ac- 
tual spectrum of the object is generally not known, because if 
it were, the required flux would be extracted from this spec- 
trum, or spectra would even be compared directly. So, we are 
dealing with some approximation of the spectrum, and we will 
now discuss how the K-correction depends on inaccuracies of 
this spectrum. 

Let us consider a supernova at z = 0.05: transforming the 
measurements from the observer filter to the same restframe 
involves a filter shift of 1 .05, while standard filters have a rela- 
tive width (width of the bandpass over central value) of about 
0.25 (see Fig 1). It means that the numerator and denominator 
of equation 6 overlap for about 80% of their values, and hence 
the correction does not depend too much on the precise knowl- 
edge of S. This is the situation that prevails for low -redshift 
samples (e.g. [121]). At z — 0.25, the overlap between T\(A) 
and Ti(A/(l + zj) becomes essentially null, but using the next 
redder filter as T 2 restores essentially a full overlap. For exam- 
ple, using the standard UBVRI filters (figure 1), W(A = 5450) 
is a very good approximation of restframe B(A = 4350) at 
z = 0.25. So, as common sense indicates, the measurement 
filter should follow as much as possible the red-shifted target 
filter. This gives rise to the cross-filter K-correction concept 
proposed in [202], where T\ and T 2 are assumed to be differ- 
ent (in the same frame). When K-corrections are expressed 
in magnitudes, transforming for different filters introduces an 
additional term accounting for the different definition of mag- 
nitudes in the involved bands. 

The worst case for K-corrections happens roughly when 
T\{A) and T 2 (A/(1 + z)) overlap by about half of their width: 
then about half of the flux in the target filter has to be de- 
rived from the assumed spectrum. In [203], a detailed study 
of K-corrections for supernovae concludes that rather than try- 
ing to obtain detailed knowledge of the SN spectrum, one 
should first pay attention to the spectrum broadband slope 
(i.e. colour), which varies from supernova to supernova. The 
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scheme proposed in [203] consists in tilting a template spec- 
trum in order to reproduce a measured colour, and then us- 
ing this tilted spectrum to compute K-corrections. One can 
note that to first order, the details of how the spectrum is tilted 
do not matter. In order to apply this tilting scheme, one has 
to measure an SN colour, which happens to be essentially 
mandatory for distance estimation. 

The Supernova Factory survey (§ 11.1) collects spectro- 
photometric SN "light curves", or photometric spectral se- 
quence, thus eliminating the need for K-corrections of any 
kind when comparing supernovae. 



8.2. Brighter-bluer relation and extinction 

Astronomical observations in the visible are essentially al- 
ways subject to some sort of extinction, generally by dust. 
For example, the Milky Way contains dust, and extragalac- 
tic observations (including supernova measurements) have to 
be corrected for Milky Way dust extinction, nowadays using 
the maps derived from dust emission in the far infrared [204]. 

The extinction by dust has been characterized mainly in the 
Milky Way [205]. By comparing colours of stars of the same 
spectral types, but affected by different line-of-sight extinc- 
tion, the authors could derive an empirical "extinction law" 
(called CCM law after Cardelli, Clayton and Mathis) which 
relates extinction at various wavelengths. The total extinction 
is found by extrapolating the found wavelength dependence 
to a null photon energy. The first obvious parameter that char- 
acterizes a line of sight will then be its optical depth, i.e. the 
amount of absorbing or scattering medium. This is conven- 
tionally parametrized (and in particular in [205]) by the red- 
dening of the B - V colour (defined in §3). Naming E(B - V) 
the colour excess attributed to dust, the extinction (in mag- 
nitude) at a given wavelength reads A(A) = E(B - V) * R(A), 
where R(A) is the "extinction law", directly related to the cross 
section of photons on the dust. Once R(A) is known, one 
may deduce the line-of-sight extinction from the reddening 
of colour. Identical objects extinguished by different amounts 
of intervening dust exhibit colours which lie along lines in 
colour-colour spaces. 

The dust extinction law should depend on the nature of the 
intervening dust and [205] proposes a one-parameter family 
of extinction curves, parametrized by Ry = R(Ay). The Milky 
Way dust is characterized on average by R v = 3.1. 

SNe la exhibit a brighter-bluer relation (see e.g. [119] and 
references therein), as one would expect from dust, and also 
exhibit a strong correlation of colour relations [109, 206], 
which constitutes a tantalizing analogy. However, this justi- 
fies parametrizing the colour variations of supernovae in the 
same way as dust extinction, but does not impose the assump- 
tion of CCM extinction laws. As we will discuss later, super- 
novae allow one to carry out a similar analysis as the one done 
to derive CCM extinction laws, in order to derive an average 
colour law for supernovae. Using then the standard candle 
property of SNe la, one can extract the equivalent of Ry from 
the supernova data. One should note that in order to measure 
an SN colour, one has to carry out measurements in at least 



two bands. 



8.3. Supernova empirical light curve models 

We are now in a position to express some minimal require- 
ments of supernova models used to "compress" the photomet- 
ric data characterizing an event: in order to derive a distance 
we need a brightness (anything that scales linearly with the 
observed flux), some measure of the decline rate and some 
measure of colour. 

If the measurements are dense enough along the light curve, 
there is no crucial need for an explicit model, which was the 
approach originally followed for the Amis paradigm [102]. 
From well-sampled light curves, one can build smooth dis- 
crete templates with different Ara^ values, where this quan- 
tity indexes the decline rate measured as the magnitude dif- 
ference between peak and 15 (restframe) days later. An inter- 
polation scheme (figure 1 of [207]) allows one to fit possibly 
sparse observations. Templates have been continuously re- 
fined [110, 208-210]. The Am^ model was in particular used 
in the Calan-Tololo survey analysis ([8, 121] , §11.1), and one 
of the acceleration discovery papers ([10], §11.2). Recently, 
the Snoopy model [211] revisited the Am\s paradigm and ex- 
tended it to NIR. 

Light curve templates for standard visible and NIR bands 
were determined in [206], mainly to assess the empirical re- 
producibility of SNe la, and could then be used to fit SN pho- 
tometry. In order to account for the decline-rate variation, the 
"stretch" paradigm [103] proposed to stretch the time axis of 
these templates in B and V standard bands (§3). As data im- 
proved in size and quality, these templates were refined and 
extended to the [/-band [104, 195]. The stretch paradigm 
works poorly for bands redder than R. The stretch approach 
was used in one of the accelerated expansion discovery pa- 
pers ([ 1 1 ], § 1 1 . 2) . All the models sketched above rely on light 
curve templates in standard bands. 

The Bayesian Adapted Template Match (BATM) method 
[212] relies on about 20 well-measured events (spanning a 
decline -rate range) and about 100 spectra to predict the ob- 
served light curves in a given instrument at a given redshift, 
for a given extinction. Using an extinction prior, it eventu- 
ally delivers a distance estimate, sensitive to the prior choice. 
It was used to extract cosmological constraints from the first 
large SN compilation [212]. 

Currently, there are two main streams of SNe la empirical 
models in the visible used to derive cosmological constraints: 
MLCS on the one hand, SALT, SALT2 and SiFTO on the 
other, and we will now describe those. 



8.3.1. The MLCS model 

The Multi Light Curve Shape (MLCS) model consists in a 
one-parameter family of light curve shapes for standard vis- 
ible bands BVRI (i.e. from 400 to 900 nm, figure 1). This 
complete model of SN light curves depends on three parame- 
ters (plus a reference date): a distance modulus, a brightness 
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offset (due to the brighter-slower relation, which indeed in- 
dexes light curve shapes), and an extinction value (in the V- 
band). At variance with essentially all other SN light curve 
models (BATM is another), MLCS incorporates the brighter- 
slower and brighter-bluer relations, which imposes to train it 
on events at known distances [95, 213]. The initial modest 
training sample of 12 events (by 1995) was later enlarged to 
about 100 events (around 2002) and the model extended to- 
wards the blue by adding the t/-band (320-400 nm) [214]. 
This second version is named MLCS2k2. Prior to fitting, mea- 
surements have to be translated into standards restframe bands 
using K-corrections. 

The MLCS(2k2) model assumes that extinctions follow 
some known distribution [95] (the "extinction prior"), which 
is not mandatory but improves the formal uncertainty of 
colour measurements (figure 4 of [95]) and enforces positive 
extinctions. The derivation of this prior in [214] is sometimes 
described as a statistical separation of intrinsic (internal to the 
supernova) and extrinsic (due to dust extinction) colour vari- 
ations, which means that there is no separation on an event 
per event basis. This statistical separation only matters for 
the shape of the prior. The other assumption of the model is 
that after accounting for colour variations correlated with light 
curve shape, the average colour relations are due to extinction 
by dust, similar to the average Milky Way dust. 

Both of these assumptions have been extensively discussed: 
first, the extinction prior can bias distances (as argued in, e.g., 
[195], appendix of [11]). Reference [200] shows from simu- 
lations that distances are unbiased if the prior faithfully repro- 
duces the actual extinction distribution, and the accuracy of 
colour measurements is accurately known. This distribution 
is, however, not known from first principles, and one should 
expect some difference between the assumed and actual distri- 
butions, but in a cosmological context, a difference only mat- 
ters if it is redshift-dependent. So, using a prior assumes at 
least that the SN population at all redshifts is drawn from the 
same parent population, possibly with a redshift-dependent 
selection function (§2.3 of [215], §5.1 of [200]). Note that 
as any other light curve model, MLCS can be used without 
any extinction prior. We will discuss the Milky Way extinc- 
tion dust hypothesis in §8.3.5. 



8.3.2. SALT and SALT2 models 

SALT [108] models supernova light curves by integrating 
a spectral template (i.e. a spectral energy distribution as a 
function of SN phase) into properly shifted observer filters. It 
implements the light curve width variability mostly using the 
"stretch" paradigm (§8.3 and [103, 104]). The training con- 
sists in reproducing photometry of a training sample by ad- 
justing a broadband phase-dependent correction to a spectral 
time series from [203]. The colour variations are modelled as 
an extinction factor exp(-cCi(A)) multiplying the flux, where 
c is the colour parameter of each supernova, and Ci(A) is the 
"colour law", common to all supernovae. It is fitted during 
the training, and found to be poorly compatible with extinc- 
tion laws. There is no assumption on the source on colour 



variations and colours that might reflect negative extinctions 
are accepted. Note that, even if colour variations induced by 
extinction are one-sided, the measured colour value is still an 
unbiased estimate of true colour (e.g. §5.2 in [216]). The 
SALT spectral template spans [300,800] nm, which covers the 
restframe [/-band, as required to fit high-redshift data from 
the Supernova Legacy Survey, for which SALT was initially 
developed. As a spectral template, SALT incorporates K- 
corrections into the model. SALT does not attempt to model 
distances, and hence any sufficiently well-measured super- 
nova can enter the training sample. SALT parameters are m* B , 
the peak brightness in restframe Z?-band, s, the stretch factor, 
and c, the (restframe) B - V colour at maximum, plus some 
time reference. 

SALT2 is a further development of SN spectral templates 
[106, 114]: the variability of spectra (and broadband light 
curves) of supernovae is modelled using principal compo- 
nent analysis (PCA), and the whole spectral model and its 
variability is fitted from a set of training spectra and light 
curves. As in SALT, the "colour law" is fitted to the data 
(and found again incompatible with known extinction laws). 
The first variability component found is extremely similar to 
stretch (or decline rate) and no significant further component 
is found. Because SALT2 does not model distances, it can in- 
corporate high-redshift events in its training sample, and do- 
ing so allows it to go as blue as 250 nm, i.e. bluer than ac- 
cessible from the ground using low-redshift events. SALT2 
delivers the average spectral model and its first variability 
component, the colour law, the uncertainties, and some es- 
timate of the intrinsic variability of supernovae around these 
averages. SALT2 finds that broadband light curve amplitudes 
scatter (coined "colour smearing" in [200]) by at least 2.5% 
around the average model, with a rapid increase when going 
bluer than ~360 nm. At variance with all other SN mod- 
els, SALT2 training actually determines spectral sequences, 
but the large number of required parameters constitutes the 
main limitation of the model. This flexibility of the training 
process imposes the necessity to resort to some regulariza- 
tion scheme in order to actually define the template in regions 
of the phase/wavelength/decline-rate space poorly covered by 
the training sample. SALT2 delivers a model uncertainty. 



8.3.3. The SiFTO model 

SiFTO [109] models SNe light curves from a SED tem- 
plate, constructed through a smooth correction of the empiri- 
cal template from [217]. The light curve shape variability is 
accounted for using the stretch paradigm, using a measured re- 
lation between stretch factors in different bands. When fitting 
a supernova event SiFTO adjusts a phase-independent mul- 
tiplicative function of wavelength in order to match exactly 
the measured light curve amplitudes. The obtained mangled 
SED is then used to evaluate the light curve amplitudes in 
the standard UBVR bands plus a bluer band, covering in total 
[250,710] nm. This set of restframe magnitudes readily pro- 
vides an overall amplitude to* , and a single B — V colour de- 
rived from an optimal combination of all measured restframe 
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colours. This combination relies on restframe colour-colour 
diagrams obtained during the training, which also delivers a 
model uncertainty, used when fitting. As SALT and SALT2, 
SiFTO does not predict distances and can then use any well- 
measured event for training. Its small number of parameters 
makes its training easy and precise. As SALT and SALT2, 
SiFTO finds colour relations (e.g. U-B versus B — V) incom- 
patible with extinction by dust, even allowing for a coloured 
variation of supernovae related to stretch. 



8.3.4. Tripp distance estimators 

SALT, SALT2 and SiFTO do not directly deliver distances. 
All propose (but do not impose) to estimate distances using 
the most general distance estimator linear in decline rate s 
(e.g. stretch) and colour c: 

fi = m* B + a(s - 1) - Be + Const (7) 

where p is a distance modulus estimate (i.e. it varies as 
5 log\o(d L ) up to an additive constant), a and B are global 
parameters that characterize, respectively, the brighter-slower 
and brighter-bluer relationships. These parameters have to be 
empirically determined, typically along with the cosmological 
parameters, as proposed by Tripp in [119]. m* B refers to the 
peak magnitude that would be observed in the filter matching 
B restframe, where the choice of B is arbitrary. For most cos- 
mological SN analyses, the measurements are actually carried 
out around the B restframe band and the distance uncertainty 
is dominated by the Be term. Assuming that colour variabil- 
ity at a fixed decline rate is due to Milky Way like dust sets 
B = R B = 4.1. 

The most common criticism to SALT, SALT2 and SiFTO is 
that the reported colour parameter mixes both intrinsic colour 
variations (slower supernovae tend to be bluer) and extrinsic 
ones (e.g. dust extinction), and the distance estimator defined 
in equation 7 conflates both effects via a single brighter-bluer 
coefficient/? (e.g. §2.1 of [215], §2 of [136]). This is subopti- 
mal if there are two or more sources of colour variations with 
different B values, but more importantly biases cosmology if 
the admixture building up the colour distribution evolves with 
redshift. Let us first assume the extra source of colour varia- 
tion is entirely related to decline rate: assuming we know the 
relation, we can then separate the contributions to colour for 
each event. To first order, this separation does not change the 
form of equation 7, but only alters the value of a (B is unaf- 
fected, see § 4.1.1 in [108], §4.2.3 in [106]). Hence, account- 
ing for decline rate related colour variations does not change 
distances evaluated using equation 7. So, this criticism as- 
sumes that there are at least two contributions to colour vari- 
ations, independent of the decline rate. This is a serious con- 
cern that can only be handled via more parameters per event 
and hence applies to all current photometric distance estima- 
tors. 



8.3.5. Comparing light curve fitters 

SALT2 and MLCS2k2 were compared thoroughly in [200], 
which finds that they deliver systematically different dis- 
tances as soon as they make use of the restframe f/-band 
([320,400] nm). This difference in distances has a dramatic 
effect on cosmological fits, and is traced to two causes: differ- 
ent predicted t/-band peak values at given B and V and the ef- 
fect of the MLCS2k2 extinction prior on high-redshift ground- 
based events from the SNLS sample (see §1 1 of [200]). When 
MLCS2k2 uses SALT2 light curve templates, and does not 
apply any extinction prior, the cosmological results become 
essentially as when using SALT2 right away. So, [200] argues 
that some unidentified issue with the low-redshift observer- 
frame [/-band seriously affects MLCS2k2, and accounts for 
the difference with SALT2 as a systematic uncertainty. On a 
similar data set, [218] identifies that observer-frame f/-band 
SN data are incompatible with higher redshift measurements 
of the same restframe spectral region, and ignores this chunk 
of data. 

Regarding restframe U-band, the difference between 
MLCS2k2 and SALT2 models lies in their training strategy: 
while MLCS2k2 has to train restframe U from observer-frame 
U SN data, SALT2 (and SiFTO) use higher redshift events 
(typically at z ~ 0.2 - 0.5) to model restframe {/-band. 
Observer-frame U is notoriously difficult to calibrate because 
the actual shape of the bandpass is determined by the variable 
atmospheric absorption 5 . This difficulty is illustrated, in the 
supernova framework, by the observed dispersion of the rest- 
frame U -band light curve amplitude, which is much larger for 
nearby events (where restframe U is measured from observer 
U) than for higher redshift events (where redshift pushes rest- 
frame U into regions unaffected by the atmosphere (see fig- 
ure 11 in [196], §2.6 in [218]). For MLCS2k2, dropping the 
observer-frame [/-band at training implies that the resulting 
model misses restframe {/-band, which is unacceptable: all 
modern SN models insist on modelling near UV SED of su- 
pernovae, because distances to many high-redshift events rely 
on this modelling. 

In [200], the second source of MLCS2k2/SALT2 difference 
is traced to the use of an extinction prior for high-redshift 
SNLS events. The fact that applying or not the extinction 
prior to a measurable quantity (the event colour) questions 
the prior (as discussed in §8.3.1), and dropping the prior is 
technically possible. We also note that for samples of several 
hundred events, the uncertainties associated with the choice of 
the prior (table 7 of [215], figure 32 of [200]) are larger than 
the statistical uncertainties. Finally, as events are better and 
better measured, the statistical benefit from constraining their 
parameters a priori is vanishing. 

Colour variations independent from decline rate are typ- 
ically assumed to be due to extinction by dust or not. In 
this second case, both colour relations (the colour law) and 



5 Correction of the absorption by the Milky Way dust might also be large in 
observer (7-band, and correspondingly uncertain [204]. 
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the brighter-bluer coefficient (6 of equation 7) should be ex- 
tracted from the data. In [200] (§11), it is shown that inserting 
the fitted SALT2 colour law into MLCS2k2 has no sizable 
effect on MLCS2k2 distances. Fitted brighter-bluer coeffi- 
cients are generally smaller than those for Milky Way dust 
(see e.g. [119] and references therein, [108]). The finding 
from MLCS2k2 distances that we might be at the centre of a 
"Hubble bubble" [214] (i.e. a spherically symmetric under- 
density) was shown to be a consequence of assuming rather 
than fitting the brighter-bluer relation [136, 219]. 

In [136], one can find a comparison of SALT, SALT2 and 
MLCS, but the SALT and SALT2 packages were probably not 
used properly: using the same photometric data, supernova 
fits differ from the ones published by the authors of the pack- 
ages. In [215] SALT and MLCS2k2 were compared and found 
to agree very well. This is not very surprising as their train- 
ing samples vastly overlap. In [106], SALT2 and SiFTO are 
compared using exactly the same event sample and the found 
differences are acceptable for a few hundred events. 

Comparing the scatter of the resulting Hubble diagram 
yields useful comparisons of distance estimation strategies. 
For large and similar samples, table 1 1 in [200] exhibits rms 
Hubble residuals of 0.21 mag (about 10% in distance) using 
MLCS2k2, while SALT2/SiFTO combined distances deliver 
about 0. 17 mag (table 4 of [218]). 

SALT2 and SiFTO are compared in [106] on exactly the 
same distant events, and SiFTO outperforms SALT2 by a 
small but significant amount, 0.150 versus 0.173 mag. 



8.4. Other distance estimators 

Most distance estimators proposed for SNe la follow the 
form of equation 7. We review here two notable exceptions. 

The so-called "Bailey ratio" [220] was optimized to obtain 
the best distance estimator from a sample of 58 photometric 
spectra from the SuperNova Factory (§11.1). It uses the flux 
ratio of two narrow spectral regions F(642 nm)/F(443 nm) as 
a single intrinsic luminosity estimator, and achieves a smaller 
Hubble diagram scatter than all other known distance estima- 
tors on the same sample. This good performance remains un- 
explained, and was confirmed using long slit spectra of nearby 
events [113]. Note that at z > 0.6, most of the current SN 
spectra do not cover the 642 nm region, although the Si line 
characteristic of SNe la is around 615 nm. 

The CMagic distance estimator [221] makes use of the post- 
maximum trajectory of SNe la in the colour-magnitude (B- 
V, B) plane, in order to evaluate the brightness when a given 
colour is reached. The obtained brightness is then corrected 
for brighter- slower relation to yield distances, which exhibit a 
reduced scatter with respect to other distances, on the nearby 
training sample. The technique was applied to a cosmologi- 
cal analysis [222] and proved to be applicable to high-redshift 
events, provided that the latter are measured in bands close to 
restframe B and V, a demanding constraint. The performance 
improvement on this sample is not large enough to compen- 
sate for the loss in event statistics due to the requirements in 
phase and wavelength coverage. 



9. PHOTOMETRIC CALIBRATION 
9.1. The traditional way 

Stellar fluxes are essentially always measured relative to 
some other star. The ratio of the studied object flux to a stan- 
dard star is commonly reported as a magnitude. In order to 
measure luminosity distances, one has to convert magnitudes 
into fluxes, relying on a star that has both magnitudes (in the 
same system as the standard) and a spectrum in physical units; 
one also needs to know the pass bands of all involved photo- 
metric measurements. 

In order to report fluxes from relative measurements, one 
obviously needs primary calibrators, in practice non-variable 
stars with known spectrum. Establishing the spectrum of 
Vega (or a-Lyrae) is described in [223], where the astronom- 
ical instrument is calibrated against artificial black bodies (or 
sources themselves calibrated to black bodies). As Vega is far 
too bright to be measured in most modern instruments 6 , an- 
choring photometry of m ^ 20 objects to Vega (m = 0) goes 
in practice (see §8.1 in [224]) through two steps: the Landolt 
star network [225, 226] (< m >- 12), referenced to Johnson 
[227] and Cousins [228] standards (< m >^ 5), themselves 
referenced to Vega [229]. Although this lengthy sequence was 
not conducted to ensure a robust conversion of magnitude to 
fluxes, all supernova cosmology works until 2006 calibrate su- 
pernova magnitudes against Landolt stars and convert those to 
fluxes using the Vega spectrum (usually from [223]), because 
they had no alternative. One should note that the Landolt sys- 
tem is not a natural system (§8.1 in [224]), meaning that the 
reported magnitudes are not directly measured, but rather de- 
rived from actual measurements and colours, and hence that 
the filters defining the reported magnitudes are unclear. The 
networks of standard stars were designed to allow astronomers 
to precisely compare magnitudes of stars, but allowing one to 
convert these magnitudes to high precision fluxes was not a 
primary goal. 

A new and shorter path to convert magnitudes into fluxes 
was designed in the framework of the HST flux calibration: 
the Calspec project 7 reports spectra measured from the HST 
and carefully cross-calibrated. The instruments used to gather 
these spectra are flux-calibrated from stellar models of pure 
hydrogen white dwarfs (references in [230]). The models 
make use of temperatures and gravity obtained from Balmer 
line shapes from ground-based high-resolution spectroscopy. 
The consistency with the calibration route sketched in the 
above section was checked by bravely measuring Vega in the 
same programme [231]. As a few Calspec standards now 
have broadband magnitudes in the widely used Landolt sys- 
tem [232], the numerous SNe la distances calibrated on this 
system now have a direct link to this spectral library. The 
overall flux scale of the Calspec library is set to match ground- 
based photometry, and dominates the uncertainty budget at 



6 A SN la at z~0. 1 typically peaks at ms — 20, i.e. 10 8 times Vega. 

7 http://www.stsci.edu/hst/observatory/cdbs/calspec.html 
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about 2% [233], but is irrelevant for cosmology with SNe la. 
Because the ratio of fluxes of supernovae at different redshifts 
depend on the ratio of fluxes of the standards at different wave- 
lengths, the "colour uncertainty" of the spectral library is, in 
contrast, extremely relevant to our subject. Two terms domi- 
nate this uncertainty: the model slope (dominated in the visi- 
ble by temperature uncertainty), and the reproducibility of the 
Calspec measurements. Both terms are estimated to typically 
0.5% or less (table 12 in [224]) for the measurement of the 
star BD +17 4708 from [234]. 



9.2. Instrumental calibration from calibrated laboratory 
sources 

Using artificial sources to calibrate instruments is not a 
modern idea (see e.g. [223] and references therein). Sug- 
gestions to calibrate imagers from artificial light sources have 
been recently formulated [235], relying on photo-diodes cal- 
ibrated by the US NIST. Results from a setup delivering a 
monochromatic illumination of a screen in the telescope dome 
are provided in [236]. However, reflections on the numer- 
ous optical surfaces in front of wide-field imagers (typically 
at least three lenses, one filter and the cold chamber window) 
contribute to the image obtained with a flat-field screen, but 
not to star photometry. For the specific case of the MegaPrime 
instrument on CFHT, direct evidence for these internal re- 
flections is provided in [237]; [224] evaluates to about 10% 
or more their contribution to twilight images. This justifies 
proposing a different artificial light source setup in [237, 238]. 
The most recent results presenting calibration from an artifi- 
cial source [239] are encouraging but do not yet outperform 
the stellar approach. 

There are great expectations from these illumination de- 
vices since most ground-based wide-field survey instruments 
either own one or plan to install one. We can sketch plausi- 
ble outcomes of current efforts roughly ordered in increasing 
difficulty and reward: daily monitoring of the instrument re- 
sponse, measuring the bandpass profiles in situ, establishing 
a flat-field correction, measuring the relative transmission of 
photometric bands, or even absolute transmission of photo- 
metric bands. 



9.3. Ground-based observations and atmospheric extinction. 

The canonical procedure to evaluate atmospheric extinction 
consists in following sources at various elevations and ex- 
trapolating above the atmosphere, or more commonly around 
zenith. One then relies on repetition to beat down atmospheric 
variability. A more precise technique, relying on spectroscopy 
and atmospheric extinction models is implemented in [240]. 

Atmospheric extinction is less critical when calibrating 
from stellar standards, and again, repeating the calibration of 
stars around supernovae over tens of nights allows one to mit- 
igate the atmospheric variability (e.g. [224]). Ground-based 
rolling supernovae searches incorporate these calibration ob- 
servations in their programme. 



9.4. Linearity issues 

Since standard stars are much brighter than tertiary stars, 
themselves often brighter than supernovae, one should wonder 
whether the reported signals scale as fluxes. One first varies 
the exposure time from typically 1 s to a few hundred sec- 
onds, in order to compare different brightnesses in the same 
electronic range. Precise shutter timing is required for CCDs 
and can be obtained (see e.g. §12.1 in [224]). The overall lin- 
earity of CCDs and their readout electronics can be checked 
on the sky (App. C in [224]), alleviating concerns that quasi- 
uniform illumination ("flat fields") might not be representa- 
tive of the instrument response to stars. CCDs in space suf- 
fer from charge transfer inefficiency that builds up with the 
integrated radiation dose, alters linearity, and should be cor- 
rected for (e.g. [195,241]). Observations with the NICMOS 
NIR instrument on HST suffer from a poorly understood non- 
linearity issue, which adversely affects fluxes, and hence dis- 
tances to high-redshift supernovae obtained with this instru- 
ment (see §5.1.3 in [218]). 



9.5. Prospects for improvements 

As calibration constitutes the leading systematic uncer- 
tainty of current supernova cosmology results [218, 242], im- 
proving over the present situation is mandatory. The weak- 
nesses of the Landolt system [218, 224] have a decreasing 
impact as new supernovae calibrated via alternative paths are 
published: for example, the SDSS [197] directly ties its mag- 
nitudes to solar analogues with CALSPEC spectra. Some pass 
bands of the CSP nearby sample [98, 243] are calibrated to the 
Smith magnitude catalogue [244], which is now anchored as 
well to CALSPEC standards. The calibration of instruments 
against laboratory standards can provide a significant break- 
through, provided that parasitic reflections in the optics are 
controlled well enough and atmospheric extinction accurately 
measured. These two complications vanish if the artificial 
light source is mounted on a satellite, as proposed in [245]. 



10. ASTROPHYSICAL SYSTEMATIC UNCERTAINTIES 
AFFECTING SN DISTANCES 

10.1. Evolution of supernovae 

Evolution of supernovae refers to the concept that super- 
novae could be different at high redshift (i.e. early in the stel- 
lar history) than at low redshift. We already know that a de- 
mographic evolution exists (§5.3), and our concern here is to 
constrain if supernovae exhibiting the same observables could 
have a different average brightness at low and high redshifts: 
the derivation of distances used for the cosmological analysis 
makes the hypothesis that it is not the case. One set of tests 
compares properties of nearby and distant samples, and the 
other set compares properties of events at similar redshifts in 
different environments, i.e. different stellar ages. 
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10.1.1. Light curve rise time 

Reference [161] proposes the approximate rule that the 
maximum (bolometric) luminosity equals the instantaneous 
energy released by f3 and y decays. This relates the rise time, 
peak luminosity and produced 56 Ni mass, so that a change 
in rise time at a given nickel mass implies a change in peak 
brightness. Rise times at a given decline rate hence constitute 
a plausible handle on intrinsic peak brightness. Even without 
theoretical guidance, the rise time of the light curve is sensi- 
tive to the explosion energy and the evolution of the opacity 
of the ejecta, and any evolution of rise time with redshift un- 
related to usual decline-rate indicators is a direct threat to rel- 
ative distances, as they are evaluated from a single light curve 
shape modelling. 

In order to measure rise times accurately, authors devise 
methods that scale the actual measurement of the rise time to 
that of the fall time, so that variations of the decline rate do not 
smear out the measurement. Three studies [89, 104, 105] do 
not detect any sizable difference between distant and nearby 
samples. One might note, however, that different studies find 
different average rise time values because of different defi- 
nitions, and/or different fiducial supernova. Reference [246] 
discusses some difficulties of the measurement and closes a 
dispute caused by the troublesome comparison of rise time de- 
terminations obtained through different estimators. Rise times 
are also compared as a function of the "stellar age" of the host 
galaxy and they turn out to be in excellent agreement (table 2 
in [89]). 



10.1.2. Spectroscopic evolution tests 

Comparisons of supernova spectra across redshifts appear 
in the very first SN cosmology papers (e.g. figure 11 in [10]). 
Quantitative comparisons involving line velocities [247], line 
velocities and profiles [248], line velocities and strengths 
[249] did not reveal compelling differences. 

More detailed comparisons with increasing statistics are 
presented in [250], where global agreement is found, with per- 
haps some small differences in line strengths, however, with- 
out subtraction of the host galaxy light. Studies of samples of 
similar size [115, 1 16, 251] do not confirm these small differ- 
ences, but they carry out galaxy light subtraction. 

High-redshift samples allow one to study the UV part of 
the spectrum (below ~350 nm) absorbed by the Earth's atmo- 
sphere for nearby events, and deemed to encode metallicity 
of the progenitor, although theoretical predictions disagree on 
the sense of the effect [252, 253]. These high-redshift stud- 
ies agree on a event-to-event scatter increasing towards blue 
wavelengths (figure 7 in [114] shows the same trend using 
photometry only), and stress that comparison at wavelengths 
below ~350 nm lack a nearby counterpart, which is nowa- 
days being collected from space. A first sample of 12 events 
[254] shows reasonable agreement in the UV part with events 
at z ~ 0.5 and mildly suggests that the observed differences 
might reflect sampling variance. The large observed variabil- 
ity in UV, probably due to metallic features [254], suggests 



that cosmological analyses aiming at higher and higher pre- 
cision limit their use of restframe UV (below ~ 350 nm) to 
measure distances to high-redshift supernovae. 

To summarize, when differences between nearby and dis- 
tance spectra are observed, they can usually be interpreted as 
the already known demographic evolution, i.e. that explosions 
tend to be more energetic as redshift grows. 



10.1.3. Segregation according to host galaxy types 

In order to compare brightnesses (or even derived dis- 
tances) of events at different stellar ages, one can take advan- 
tage that varying stellar ages cohabit at the same cosmolog- 
ical age. For example, the local Universe hosts both young 
stellar systems (and even galaxies producing stars), as well as 
old evolved heavy elliptical galaxies. One can then question 
whether we derive identical average distances to events at the 
same redshift but in different environments, thus directly ad- 
dressing the heart of the evolution question without strong as- 
sumptions on the distance-redshift relationship. Comparisons 
between spiral and elliptical galaxies reported null differences 
(e.g. §5.3 in [199], [255]). Many works test the effects of host 
galaxy metallicity (see §1. in [256]), and report null results, 
with the exception of [257] which only tests a small sample 
of early-type hosts. A first indication of systematic difference 
of distances as a function of the host galaxy stellar mass was 
proposed in [258], and confirmed on independent samples by 
[128, 132], and could perhaps be attributed to systematic dif- 
ferences in metallicity. Note that the galaxy stellar mass is 
inferred from the wide-band photometry of the host galaxy, a 
by-product of light curve measurements, and hence the effect 
is easily corrected for. Ignoring this brightness difference of 
~8% between low and high stellar-mass hosts would bias the 
dark energy equation of state w by ~0.04 to 0.08 (table 9 in 
[132]). 



10.2. Gravitational lensing of supernovae 

Gravitation bends light rays along the transverse compo- 
nent of the gravitational potential gradient (e.g. [259, 260]). 
Mass concentrations act somehow as lenses, which coined this 
class of phenomena. Cosmic bending of light rays re-maps the 
image plane, which alters the brightness, shape, size and posi- 
tion of sources. In the framework of supernova observations, 
the gravitational lensing of supernova light by foreground 
mass contrasts induces two main effects: an increased scat- 
ter of the observed brightness, and when going through strong 
lenses, the possibility of multiple images of the same event, 
usually with different time delays. Strong lenses are character- 
ized by a singular mapping of the source plane onto the image 
plane, leading to an infinite magnification on critical lines, and 
the possibility of multiple images (see e.g. [259, 261]). The 
lines of sight subject to critical lines are very scarce and the 
vast majority of supernovae incur asymptotically weak lens- 
ing. Because supernovae are point sources, their geometrical 
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shape is unaffected by lensing, but their brightness follows the 
optical magnification. 

Lensing of supernovae causes an extra scatter of their ap- 
parent brightness: the rare events observed through a mass 
concentration will appear brighter, while the bulk of the events 
(going through under-dense regions) will be slightly fainter 
than the average. Lensing conserves the average intensity 
[262], and does not constitute a first-order threat on lumi- 
nosity distance measurements [263]: the mass density that is 
probed through these measurements is the actual average den- 
sity, rather than that of a smooth component 8 . Strong lensing 
might violate flux conservation in the case of several images 
where one or more escapes detection. This is commonly in- 
vestigated through numerical simulations (e.g. [264-266]), 
and as expected, the result strongly depends on the assump- 
tions on the "dumpiness" of matter distribution. Monte-Carlo 
techniques also face practical difficulties in re-summing mul- 
tiple images that would be merged by an observer [265, 266]. 
Fortunately, one may use systematical searches for multiple 
images of distant radio sources [267], in order to evaluate an 
upper bound to the average amount of flux lost in faint sec- 
ondary images of supernovae: in [196] (§7.6), it is found that 
at z — 1, less than 0.3% of the flux is lost to secondary im- 
ages. This is far below other systematic uncertainties, both 
for present and foreseen SN surveys. Note that for very large 
statistics SN surveys, the frequency of multiple images be- 
comes measurable and this very small effect could even be 
estimated and corrected for. 

The expected magnification distribution is significantly 
skewed [262, 264-266, 268, 269], because the magnification 
decays very rapidly with the impact parameter to a lens. So, 
even if the average flux is not altered by lensing, the aver- 
age magnitude is, and its distribution is no longer Gaussian: 
least-squares estimates from high-statistics Hubble diagram 
fits made in magnitude space can be biased. In [269], it is 
shown that the maximum likelihood (rather than least squares) 
delivers unbiased cosmological parameters, even for several 
thousands of events up to z ~ 1.7. 

The extra brightness scatter due to lensing also broadens 
distance residuals and degrades the cosmological inferences. 
Again, different simulations using different assumptions on 
the dumpiness of matter lead to different conclusions: while 
[270] find an extra magnitude scatter of ~ 0.09 x z, [265] find 
0.04 at z — 1 for smooth halo profiles. 

Evidence of weak lensing of supernovae by foreground 
structures was obtained recently, by correlating the Hubble 
line residuals of individual events with the magnification ex- 
pected from the foreground galaxies [271], and also shows 
that the amplifications agree with the expectations from mea- 
sured mass contrasts. Stronger detections are only expected 
from increased statistics. From the same sample of events, 
[272] extracts weak constraints on dark matter halos. We also 
have evidence for a supernova at z — 0.6 experiencing strong 



The most common parametrizations of gravitational magnification are ob- 
tained perturbatively, thus miss the large-angle de-magnification, and hence 
do not conserve the average brightness. 



magnification by ~3. 6 (1.4 mag.) through a foreground galaxy 
cluster at z = 0.18 [273]. 

Lensing of supernovae will not probably limit the cosmo- 
logical accuracy of future high-statistics and high-precision 
SN surveys. For large-scale SN surveys (typically target- 
ing 10 4 to 10 5 events) lensing will open two new avenues: 
constraining the shape of dark matter halo from correlation 
of Hubble residuals with properties of the line of sight (see 
[271, 272] and references therein); large surveys will also de- 
tect multiple images of strongly lensed events, where an ac- 
curate measurement of the distance, flux ratio and time de- 
lay of the images provides valuable cosmological information, 
assuming the lenses can be accurately modelled (e.g [274] 
and references therein). One should note here that any kind 
of high-redshift time-variable source is eligible for this type 
of study, typically any type of bright enough supernova, and 
to a lesser extent quasars. However, the "standard candle" 
property of SNe la allows one to determine the lensing mag- 
nification factors (rather than their sole ratio), which breaks 
some degeneracy between the lens potential and the Hub- 
ble constant [275]. Simulations conducted in [276] indicate 
that the prospect for finding multiply imaged supernovae are 
modest for a possible space mission (about 10 events), where 
the ground-based LSST project would find (9(100) events in 
its 10-year mission; the same space mission would measure 
time delays of about 20 multiply imaged quasars, while LSST 
would measure several thousands. 



11. SN SURVEYS AND THEIR COSMOLOGICAL 
FINDINGS 

We review here the main SNe la surveys. Table I lists the 
main event samples we will present in the following sections. 
We will start with nearby SN samples (typically z < 0.1) 
used in essentially all analyses. Most analyses presenting new 
events also incorporate events from pre-existing samples. For 
exhaustive lists of SNe la samples, we refer to compilations, 
rapidly presented in § 1 1 .6. 



11.1. Nearby samples 

"Nearby" here broadly refers to samples at z 0.1, but 
at z > 0.01, because smaller redshifts are too contaminated by 
Doppler shifts due to proper motions. Practically, the photom- 
etry and spectroscopy of these events are mostly carried out 
using respectively, 1 and 2 m class telescopes. There are no 
practical examples of cosmological constraints derived with- 
out a nearby sample, because it increases the redshift lever 
arm, and decouples the nuisance parameter £H% (see §2.3) 
from cosmological parameters. 

While wide-field imagers make it easier and easier to find 
and measure distant supernovae, nearby events are still diffi- 
cult to find. As discussed above, the Calan-Tololo nearby sur- 
vey delivered the first well-measured sample [9] found from 
photographic plates. The Lick Observatory Supernova Search 
(LOSS) [283] has been a key provider of nearby events for 
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Name 


z range N evmts 


Dates N ba „ ds 


§ References 


Calan-Tololo 


<0.1 29 




90-93 


4 


11.1 [8, 9, 121] 


HZT 


[0.2,0.6] 10 




95-97 


2 


11.2 [10, 199,277] 


SCP 


[0.2,0.8] 42 




92-97 


1-2 


11.2 [11] 


SCP (HST) 


[0.4,0.8] 11 




97-00 


2 


11.3 [195] 


SN-Goods (1&2) 


[0.3,1.5] 37 




02-05 


2-3 


11.4 [29, 278] 


SNLS (1&3) 


[0.2,1.0] 242 ( 


-400) 03-08 


3-4 


11.5 [106, 196, 218, 242] 


ESSENCE 


[0.3,0.7] 60 (~ 


200) 


02-06 


2 


11.5 [215] 


SN-SDSSII 


[0.1,0.4] 103 ( 


-500) 05-07 


3-5 


11.5 [200,279] 


CFA(1,2&3) 


<0.1 185 




01-08 


4-5 


11.1 [280-282] 


CSP nearby 


<0.1 85 (129) 


04-09 


5-8 


11.1 [98,243] 



TABLE I. SN surveys presented in the following sections. N evenls is the number of reported events, possibly followed by the anticipated size of 
the final sample, in parentheses. Ntands indicates the number of different bands these events were measured in. A fair fraction of CSP events 
was measured in NIR, which makes 8 bands in total. 



about 15 years, and mostly targets galaxies not in the Hubble 
flow, but also delivers events in the Hubble flow. 

The Harvard Center for Astrophysics (CfA) has produced 
three releases [280-282] of spectra and light curves of nearby 
events discovered by both professionals and amateurs. With 
185 events, the last release (CfA3) constitutes the largest ho- 
mogeneous sample to date, with systematic calibration uncer- 
tainties evaluated to 3% (a probably pessimistic value accord- 
ing to [218]). The CfA data set has recently been comple- 
mented by a massive release of spectra [284], which will sig- 
nificantly contribute to training of light curve fitters. 

The Carnegie Supernova Project (CSP) also follows up de- 
tections from professionals and amateurs. It released two 
batches [98, 243] of high-quality light curves reaching a to- 
tal of 85 events among which 70 have NIR light curves. 

Regarding distance measurements the merits of the various 
samples are evaluated in detail in [218] (§2.3 and 5.1.2), and 
turn out to be comparable after strict quality cuts. The issues 
related to the U-band of nearby supernovae are discussed as 
well in the same work (§2.6) and the conclusion is that the U- 
band measurements of the CfA and the CSP are incompatible, 
and the latter agree with colours of higher redshift events. 

The SuperNova Factory aims at collecting several hundreds 
of spectro-photometric spectral sequences of SNe la [285], 
using its own blank-field search. The reduction of spectra is 
very involved and sequences have not yet been released. A 
low-scatter distance estimator relying on photometric spectra 
has been proposed (§8.4 and [220]). 

When events are used for cosmology, one should evaluate 
the distance bias caused by the search and identification steps 
possibly missing faint events (the so-called Malmquist bias). 
For a nearby sample, only the average distance bias matters 
for cosmology, and when caused by a sharp flux limit, its value 
does not depend on the value of this limit and only mildly 
depends on the shape of the cut (§7.4 in [196] and §2.7.1 
in [218]). Different works find slightly different bias values, 
from 4% ([11]) down to 2.7% ([196])) mainly because of dif- 
ferent assumed values of the brighter-bluer relation. Whether 
this flux limit model applies to a given sample can be tested 
through the phase of the first measurement of the found events 



[11,218]. 



11.2. The SCP and HZT pioneers: the twin papers 

In the mid 1990s, two groups embarked on the measure- 
ment of distances to high-redshift supernovae (typically z ~ 
0.5), the Supernova Cosmology Project (SCP) and the High-Z 
team (HZT). The first campaigns allowed the teams to refine 
the search strategy, and the "batch discoveries" of 1995 (§7.1) 
define the start of systematic programmes. From there on, 
time allocation committees awarded observing time to both 
teams on a regular basis. Both teams even gained access to 
the HST for photometric follow-up of their supernovae. 

The SCP produced an early analysis of 7 high-redshift 
events [103], which when compared with the Calan-Tololo 
nearby sample mildly favoured a matter-dominated Universe 
assumed to be flat. The quality of distances was very poor, 
and some events lacked spectroscopic identification. The HZT 
presented an early sample of 4 distant events among which 3 
are measured by the HST [286], and found £l M - 0.3 ± 0.3. 

Both teams were in fact building up larger samples, and 
published their twin analyses [10, 1 1] a few months apart. The 
SCP analysis [11] presents 42 distant events, mostly with poor 
colour measurements, and uses 17 nearby events (from the 
Calan-Tololo survey [9]). The used distance estimator ignores 
colour, however, after eliminating reddened events and having 
checked that the distant and nearby samples have the same av- 
erage colour and understood spreads. The HZT analysis [10] 
makes use of 10 distant events and has colour measurements 
for 9 of them, but colours are used through an extinction prior, 
which limits their impact 9 . It also makes use of 27 nearby 



9 The average distance modulus uncertainty without a prior would be ~0.38 
from the uncertainty of peak magnitudes, while the average reported the 
uncertainty of distance modulus is 0.18 (without intrinsic dispersion). This 
together with the fact that more than half of the events hit the prior hard 
limit Ay = indicates that the weight of the colour measurements in the 
distance estimation is modest. 
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events mostly from the Calan-Tololo survey, and finds similar 
results using two different distance estimators (which, how- 
ever, share the same extinction prior). Both distant samples 
exhibit similar photometric quality in restframe B-band with 
an average photometric uncertainty of the light curve ampli- 
tude of 0.07 mag. The SCP sample has poor or even missing 
restframe V-band measurements. Those from the HZT com- 
pare well with their restframe B-band, but however do not de- 
liver a precise enough colour that might enter a distance esti- 
mator without a prior. A Hubble diagram combining both sets 
is shown in figure 6. 




0.01 0.10 1.00 



z 

FIG. 6. Combined Hubble diagram (logarithm of a distance versus 
z) from the twin papers (top), and residuals to a Q. M = 0.3, C2a = 
cosmology (bottom). The fit favours the thick line over the two 
other proposed hypotheses. Note that the HZT events have more 
precise distances than those from the SCP. Reproduced from [287], 
with permission. 

Both teams fitted their data to an £l M , Q A Universe and 
found extremely similar best fits and contours: Prob(Q.A > 
0) > 99%, a flat matter-dominated Universe was excluded at 
more than 7 cr, Q.m — 0.28 + 0.10 for a flat Universe. Fitting 
a matter density with Qa = drives Q.m negative for both 
analyses, a certainly uncomfortable situation, which strongly 
suggests acceleration. Both teams reported a long study of 
systematic uncertainties, and did not find room for large off- 
sets: the SCP reported a systematic uncertainty roughly half 
of the statistical one. 

The significance of a non-zero cosmological constant was 



in fact modest, typically 3 cr or less, but there were two of 
them, however statistically correlated because they share a 
great part of their nearby samples. The significance of ac- 
celeration is even smaller. But taking the confidence levels at 
face value, matter-dominated universes with Q.m > 0.2 were 
excluded at more than 4 cr by each analysis. 

Both teams also explored variants of the cosmological con- 
stant by replacing it by a fluid X with unknown equation of 
state w (i.e. w is the ratio of pressure to density, see §2.1), so 
that px(z) = px(0)(l + z) 3(, " +1) . In a flat Universe, both anal- 
yses [11, 277] found similar and almost degenerate contours, 
which can be roughly approximated by £lxw = (1 - Qm)w — 
-0.72 ±0.1. 

These twin papers are usually regarded as the original key 
evidence for an accelerated expansion, and Saul Perlmutter 
(SCP), Adam Riess and Brian Schmidt (HZT) were awarded 
the 201 1 Nobel Price in Physics for this discovery. 



11.3. First Hubble Space Telescope samples 

Both pioneering analyses illustrate that the colour measure- 
ment quality is the Achilles' heel of distance measurements, 
because both teams could not use the bare colour measure- 
ments in their distance estimation. Measuring accurate light 
curves and accurate colours with the HST seems an efficient 
response to this limitation and HST photometric follow-up of 
ground-based searches started in 1997. One event from the 
SCP sample has HST measurements and three events from 
the HZT sample [286]. 

The SCP team delivered [195] the first cosmological analy- 
sis in 2003 where distances to high-redshift events use bare 
measured colours, from a mixture of 11 new HST well- 
measured events and ground-based ones. Interestingly, three 
ways of inserting colours in the distance estimator are com- 
pared (figure 9 of [195]): 

1. ignoring colours on a low reddening sample (as in the 
SCP original paper); 

2. applying a one-sided extinction prior (as in the HZT 
original paper); 

3. using bare measured colours. 

One realizes that with the third approach, the events with 
HST colour measurements totally dominate the obtained cos- 
mological constraints, which are slightly worse than the ones 
drawn in each of the twin papers. The "ground-based" events 
alone deliver very poor constraints if their bare colour mea- 
surements are used. 

Among these different approaches, one is not necessarily 
right and the others wrong; the different obtained cosmologi- 
cal constraints rather reflect different hypotheses. In the first 
method, one should check (as done in SCP original analy- 
sis [11]) that the average colour of distant and nearby events 
agrees, but if the uncertainty in the difference between the 
average colour of both samples is propagated to the end of 
the analysis, one gets the same uncertainties as with the third 
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method. So, not propagating this uncertainty relies on the as- 
sumption that both samples have the same average colour, and 
that the measured difference reflects sampling variance, pro- 
vided however that it is compatible with zero. The second 
method assumes that the colour distribution itself does not 
vary with redshift, i.e. that events at all redshifts are drawn 
from the same parent population (possibly with selection bi- 
ases, see §8.3.1). So, the first method assumes identical colour 
averages, the second identical (or known) colour distributions 
and the third makes no assumption. The third method is ro- 
bust to unknown redshift-dependent selection biases or de- 
mographic evolution, which should generally be expected to 
some level. 

This 11 -event sample of [195] also illustrates wavelength 
coverage issues: these events were measured in R and / filters 
(namely F675W and F814W of WFPC2), which map almost 
exactly on restframe B and V at z — 0.5 but on U and B at 
z ~ 0.75. The colours of three events at z > 0.7 could only 
be measured by a comparison with U - B of nearby events, 
but there were very few nearby events to compare, and the 
restframe U - B versus B - V colour relation was logically the 
dominant uncertainty of the study (table 9. of [195]). Later, 
the U - B versus B - V restframe relation was first established 
in [108] from nearby events from [281], and then using distant 
events [109, 114, 196]. 



11.4. HST searches 

In 2002, the Advanced Camera for Surveys (ACS) was in- 
stalled on HST, and opened the possibility of finding faint 
high-redshift events with HST. The first harvest of supernovae 
detected with this instrument reported 16 events [29] (comple- 
mented by 2 from [288]) among which 8 at z>l. The identi- 
fication came from a mixture of ground-based spectroscopy, 
ACS low-resolution slit-less spectroscopy, and typing from 
light curves. The highest redshift light curves are measured 
in the NIR using NICMOS aboard HST, using the F110W 
([0.8, 1 A\ym) and the F160W ([1 .4, 1 .S]fim) filters. This work 
gathers all SN la events available by late 2003, selects a "gold" 
sample and estimates distances in a consistent way, using 
MLCS2k2 (§8.3.1) with an extinction prior. Taking advan- 
tage of the increase in the redshift lever arm offered by the 
z > 1 events, the key conclusion of the paper is to provide 
evidence for past deceleration, a distinctive feature of dark 
energies similar to A. The analysis thus eliminated some dust 
models as alternatives to the cosmological constant. A stan- 
dard cosmological fit yield an the equation of state parameter 
w very close to -1 (i.e. A) with a precision of ~ 0.16. 

A second round with 22 more events was published 3 years 
later [278], using a better characterization of NICMOS (the 
key instrument for events at z > 1), establishing that the pho- 
tometric response of the instrument is non-linear. Uncertain- 
ties on corrections are not propagated in the analysis, but how- 
ever cause sizable distance uncertainties common to all high- 
redshift events (see e.g. §5.1.3 of [218] for details). The cos- 
mological analysis concentrates on variable equation of state 
fits, concludes that the data do not require w + — 1, and also 



indicates that w(z > 1) < at ~98% CL. However, the deter- 
mination of w(z > 1) still has a variance of order 1 or more. 

More HST events have been collected or published in the 
last few years within different programmes. The six events 
(two at z > 1) presented in [289] have HST data. The HST 
Cluster Supernova survey targets high-redshift galaxy clusters 
in rolling search mode, in order to increase the SN yield with 
respect to blank fields [290]. Recently, the programme deliv- 
ered accurate distances to 14 SNe [291] among which 10 are 
at z > 1. The CLASH+CANDELS SN programme runs a 
rolling search with the new WFC3-IR imager, targeting high- 
redshift events mostly in galaxy clusters. It recently delivered 
an accurate distance to an event at z = 1.55 [292]. 

Although the redshift coverage of HST samples is unique, 
they now have a modest weight in cosmological fits of the 
latest compilations, which typically gather several hundred 
events (see e.g. [242], and §11.6). 



11.5. Ground-based rolling searches: ESSENCE, SNLS, and 
SDSS-II 

The ESSENCE, SNLS and SDSS-II SN surveys are second- 
generation "high-redshift" surveys. They all aim at a signif- 
icant increase in both the amount and the quality of well- 
measured high-redshift SNe la, in order to constrain the dark 
energy equation of state. They all worked in the rolling 
search mode (§7.1), i.e. image repeatedly the same point- 
ings and extract detections and photometry of the events from 
the same image sequence. All three surveys benefited from 
large guaranteed multi-year allocations for their imaging data. 
They spectroscopically identified their candidates and mea- 
sured redshifts using 4 and 8 m class telescopes. 

The ESSENCE survey [293] ran at the CTIO 4-m (Northern 
Chile) on the 0.36 deg 2 Mosaic-II imager. It monitored in R 
and / bands 36 pointings (i.e. ~ lOdeg 2 ) every fourth night, 
for a 3-month observing season, over 5 years (2003-2008). 
The 2002 season was essentially lost to weather and technical 
problems. 

The SNLS survey [196] ran at the CFHT 3.6 m in Hawaii, 
equipped with the Megacam 1 deg 2 imager. It monitored 4 
pointings in g,r,i and z bands every fourth to fifth night, as 
long as they remained visible from 2003 to 2008. The obser- 
vations were part of the CFHT Legacy Survey 10 (CFHTLS) 
and were aiming also at galaxy studies in the very deep 
stacked images. 

The SDSS-II Supernova Survey [197] used the SDSS 2.5- 
m telescope at Apache Point (new Mexico), with its original 
1.52 deg 2 imaging camera. It monitored a 300 deg 2 equatorial 
stripe every second night for 3 months per year in u,g,r,i and 
z bands. 

The deepest survey is SNLS, which delivered usable light 
curves at z > 1, thanks to ~40 mn exposures in the red bands, 
with good IQ. Then comes ESSENCE, which used 200- and 



http://www.cfht.hawaii.edu/CFHTLS 
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400-s exposures at a similar cadence, and is thus limited to 
z H 0.7. SDSS took advantage of its fast scanning camera 
and used 54 s exposures and thus efficiently addresses the 
0.1 < z < 0.4 "redshift desert" where regular point-and-stare 
imaging with large telescopes is inefficient. ESSENCE and 
SNLS use similar instrumentation and have made opposite 
choices: ESSENCE went for a wide survey in order to obtain 
the smallest statistical uncertainty on w [293], while SNLS 
was limited to four deep pointings by the observing time allo- 
cation mechanism. It turned out that both surveys found about 
twice as many events than they could possibly identify within 
their sizable spectroscopic observing time allocations. 

SNLS has produced its first [196] and third [106, 218, 224, 
242] year analyses. The first year analysis gathers 44 high 
quality nearby events and compares those to 71 SNLS distant 
events, using distances estimated with SALT (§8.3.2). Con- 
straining supernova restframe U - B versus B — V colour rela- 
tion is performed using the distant sample (measured in more 
than 2 bands), where it is found to be much better defined 
than in nearby samples. A constant equation of state fit to 
this data set yields w = -1.02 ± 0.09(stat) ± 0.054(syst), a 
result compatible with A and systematic uncertainties domi- 
nated by photometric calibration issues. The third year SNLS 
analysis uses a much larger sample of strictly selected events 
(table 1. in [218]): 123 nearby events, 242 SNLS events, 93 
from the first year SDSS sample, and 14 HST events (see 
figure 7). The systematic uncertainties correlate events and 
the introduced covariances are propagated to the cosmolog- 
ical fits, and delivered with the data [218]. The light curve 
fitters (SALT2 and SiFTO §8.3.2 & 8.3.3) training uses dis- 
tant events in order to model the UV behaviour of supernovae 
[109, 114], which improves the distance scatter at z > 0.8 
by about a factor of 2. For a flat Universe, a constant equa- 
tion of state fit to the supernovae data, CMB, and BAO yields 
[242] w = -1.068 ± 0.08(staf + sys), where statistics con- 
tributes 0.054. Again, systematics are totally dominated by 
several photometric calibration limitations (table 3 of [242]). 
Some could be overcome in a joint analysis with the SDSS 
full sample. Fits allowing a variable equation of state w(z) = 
wq + zw a /(l + z) yield w„ = -0.98 ±1.1 compatible with 0. 
These cosmological fits are fully compatible with the cosmo- 
logical constant hypothesis (figures 8 and 9). 

ESSENCE has produced its fourth year analysis [215, 293] 
which presents 60 new events at z < 0.7. These events are 
analysed using MLCS2k2 (§8.3.1), and considerable care is 
devoted to choosing the prior 11 , which incorporates provision 
to compensate selection biases, because the images are shal- 
low. A constant equation of state fit yields w = -1.05 ± 
0.12(ifaf) + 0.l3(syst). The uncertain choice of the extinction 
prior dominates systematics (0.08), followed by the possible 
existence of a "Hubble bubble" (0.06) and the restframe U — B 
versus B — V relation (0.06). These three main contributions 
could be reduced. 

The SDSS-II SN survey has produced its first year analysis 
[197, 200] with 103 new events at 0.04 < z < 0.42. The pre- 



The same data set is analysed with a different extinction prior in [200] 



sented analysis gathers nearby samples, the SNLS first year, 
ESSENCE and HST samples. As already discussed (§8.3.5), 
the core of the cosmology paper [200] is a comparison of the 
outcome of MLCS2k2 and SALT2 frameworks, which finds 
alarmingly different cosmologies. We just repeat here (§8.3.5) 
that the main causes of the discrepancy are identified as the 
restframe U - B versus B - V relation (likely induced by pho- 
tometric calibration issues of the training sample), and the ef- 
fect of applying or not an extinction prior to the SNLS sam- 
ple. Both can be settled (and indeed are), and hence do not 
threaten the use of supernovae as distance indicators. This de- 
tailed and alarming comparison of approaches was sometimes 
interpreted as a general failure of supernova cosmology. Al- 
ternatively, one might argue that both the ESSENCE [215] and 
SDSS [200] analyses have convincingly demonstrated that un- 
certainties associated with the choice of an extinction prior are 
now too large for the present sample sizes. 

Distances to 35 SDSS, ESSENCE and SNLS supernovae at 
0.2 < z < 0.7 were measured in the restframe I-band (i.e. in 
NIR for observers) and with different systematics, confirm the 
findings of the parent surveys [216]: w = -1.05 + 0.1 3>{stat) + 
QW{sys). 

None of these second-generation surveys has finalised its 
full-sample analysis yet, but efforts on cross-calibrating the 
surveys are on-going, in order to tighten the relative flux 
scales. Because of their multi-band coverage, SNLS and 
SDSS are expected to produce a precisely cross-calibrated 
sample. In this scenario, nearby samples will not benefit from 
an improved SNLS-SDSS cross-calibration and their weight 
in the determination of cosmological parameters will be di- 
minished. 



11.6. Supernova compilations 

There is no example of a single survey delivering interest- 
ing supernova constraints on its own: covering a large red- 
shift range requires different instruments generally operated 
within different projects. Compilations of supernovae pri- 
marily target the production of cosmological constraints and 
mainly consist in producing a coherent set of distances, pos- 
sibly with their systematic uncertainties. In most of the cases, 
compilations present some new events. By delivering state-of- 
the-art and homogenized collections of events, these compila- 
tions receive a lot of attention from users primarily interested 
in cosmology constraints. 

In [212], a new distance estimation method (BATM) is 
compared to MLCS2k2 and two versions of Am^, and found 
to be "generally consistent". Using a combined distance, 172 
events are then placed on a Hubble diagram, and the cosmo- 
logical conclusions strengthen earlier findings. In [29], events 
found and measured with the HST (see §11.4) are merged 
into a large compilation, ranking events into gold and silver 
samples, analysed using MLCS2k2. The Union sample [294] 
consists of 307 events selected from more than 10 sources and 
analysed using SALT. The analysis accounts for systematic 
uncertainties, in particular photometric calibration, and deliv- 
ers correlated uncertainties. This compilation was augmented 
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FIG. 7. One of the latest SNe la Hubble diagrams (logarithm of distance versus z), which gathers well-measured events by 2010. The lower 
plot displays residuals to the best fit cosmology. A comparison with figure 6 shows that distance uncertainties and scatter of high-redshift 
events were almost halved in about 10 years. Reproduced from [218], with permission. 



into the Union 2 compilation [289], mainly adding the SDSS 
(§11.5) and CFA3 (§11.1) samples, and reaches 557 events 
fitted using SALT2, again evaluating systematic uncertainties 
in a correlated fashion. The Union 2.1 compilation [291] adds 
23 events, in particular 14 HST events among which 10 are at 
z > 1. The Union (1) compilation was also augmented into 
the Constitution compilation [136] by the addition of the CFA 
3 nearby SNe sample, analysed using MLCS2k2, SALT and 
SALT2. As noted in §8.3.5, the comparison of light curve fit- 
ters conducted there is very likely affected by an improper use 
of SALT and SALT2. 



12. PROSPECTS FOR SUPERNOVA COSMOLOGY 

The fact that statistical and systematic uncertainties con- 
tribute similar amounts to current cosmological constraints 
from supernovae forces new supernova survey projects to rely 
on new setups. Three fronts are usually simultaneously pur- 
sued: increasing the number and quality of events, improv- 
ing the quality of the photometric calibration and collecting 
higher redshift events. 

One should note that essentially all supernovae events in- 



volved in cosmological analyses to date have been spectro- 
scopically confirmed. Given that samples already on disk 
amount today to about 1000 events usable for cosmology, fu- 
ture SN imaging surveys should target at least a few thousand 
well-measured events. Spectroscopic identification of such a 
large sample lies beyond likely allocations on large-aperture 
telescopes. Recent practical studies of photometric identifi- 
cation (i.e. from light curves) propose techniques to reach 
a purity of about 95%, with an efficiency still greater than 
that of spectroscopic identification [295-297]. The required 
quality of redshifts is still debated, and "after the fact" multi- 
object spectroscopy of host galaxies constitutes the most ob- 
vious way to secure high-quality redshifts. So, broadband 
photometry of supernovae should aim in the future at provid- 
ing both distances and identification, and possibly even red- 
shifts [298, 299]. Regarding dark energy constraints, all cur- 
rent wide-field imaging large-scale projects aim at combin- 
ing probes, typically cluster counts, weak shear correlations, 
BAOs, and distances to supernovae. Two reports [300, 301] 
have reviewed in detail the scientific and technical merits 
of these probes and note that distances to supernovae is the 
most mature. Both reports also stress that the potential reach 
of each probe at constraining dark energy primarily depends 
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FIG. 8. Statistical constraints from supernovae alone (figure 7) on 
(Qm,£2a) cosmological models, including systematic uncertainties. 
Reproduced from [218], with permission. 
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FIG. 9. Statistical joint constraints on (Q.m, Vf) flat cosmological mod- 
els from supernovae of figure 7 (blue contours including systematics), 
CMB and BAOs (green). The joint contours (grey) favour models 
close to a cosmological constant, i.e. w = -1. Reproduced from 
[242], with permission. 



on systematic uncertainties, which especially for forecasts 
are highly uncertain, and strongly suggest a multi-probe ap- 
proach. Coming back to results, one might note that five years 
after these reports, dark energy constraints are still dominated 
by supernova distances. 



12.1. Space missions concepts 

Although several mission concepts to measure supernovae 
from space were developed in the last decade, none has yet 
passed the selection processes. Working from space is manda- 
tory at z > 1, because reliable distances should then be mea- 
sured in the NIR, in order to allow a direct comparison with 
nearby events measured in blue bands. Ground-based mea- 
surements in the NIR suffer from both a very bright atmo- 
spheric emission and a sizable and time-variable atmospheric 
extinction. Space also offers a very good image quality and a 
stable environment, which both significantly improve photo- 
metric measurements with respect to ground. 

Euclid at ESA is an approved dark energy project, to be 
launched in 2019, but has no compelling supernova pro- 
gramme yet: regarding dark energy, the project concentrates 
on measuring galaxy shapes for weak lensing and conduct- 
ing a (spectroscopic) galaxy redshift survey in NIR for BAOs 
[47], on about one third of the sky. It will also deliver 
NIR imaging in three bands of the surveyed fields. The 
WFIRST project at NASA is not yet approved and contains a 
high-redshift supernova programme, targeting a few thousand 
events. 



Project 


Mirror 


Area First Larj 


;e survey 




(m) (deg 2 ) light 


(nights) 


CFHT/Megacam 


3.6 


1. 2002 


500 


Pan-STARRS 


1.8 


7. 2009 


>1000 


Blanco/DEC 


4.0 


3. 2012 


500 


Subaru/HSC 


8.2 


1.8 2013 


-500 


LSST 


6.5 


10. 2019 


3500 



TABLE II. Key figures for the major past and future wide-field imag- 
ing facilities. The SNLS survey was part of the CFHTLS large sur- 
vey. Pan-STARRS is currently constructing a second telescope and 
aims at eventually operating four of them. The 500-night survey to be 
run at the CTIO-Blanco is called Dark Energy Survey (DES). LSST 
real diameter is 8.2 m but suffers from a 5 m central occultation. 
This facility will almost entirely observe in survey mode during its 
anticipated 10-year lifetime. 



12.2. Ground-based SN projects 

Planned ground-based SN surveys will enlarge the redshift 
span of Hubble diagrams (in particular compared with SNLS) 
by using thick high-resistivity CCD sensors (e.g. [302]), 
which typically deliver quantum efficiencies above 50% at 1 
//m (while the silicon cutoff lies at 1.1 /mi). Ground-based 
projects are including the four canonical dark energy probes 
in their science agenda. High-statistics Hubble diagrams with 
a high efficiency up to high redshifts (z ~ 1) is reserved to 4 
and 8 m class telescopes. Table II summarizes the key figures 
of these projects. 

For DES, the promises of the SN survey were published re- 
cently [303]: the available observing time will allow one to 
measure distances to about 4000 events at z < 1 .2, and ten- 
tatively relying on photometric redshift of both host galaxies 
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and supernovae to build a Hubble diagram. 



13. SUMMARY AND CONCLUSIONS 

During the 1990s, distances to SNe la allowed observers to 
map the expansion history of the Universe, and led to the dis- 
covery that this expansion is currently accelerating. Dark en- 
ergy has been inserted in the standard cosmological model as 
the fluid causing the acceleration, which has been confirmed 
independently of supernovae. Second-generation supernova 
surveys have delivered precise constraints on the "nature" 
of dark energy: the cosmological constant paradigm, corre- 
sponding to a dark energy of constant density, describes the 
data very well and leaves little room for dark energy density 
evolution. The ACDM cosmological model successfully de- 
scribes large-scale cosmological observables, including those 
involving structure formation. 

The nature of SNe la explosion is still questioned, and both 
progenitor scenarios have supporting evidence and face diffi- 
culties. The paradigm of a WD exploding when reaching the 
Chandrasekhar mass fulfils observational constraints, but is 
accepted by default of an alternative rather than observational 
evidence, the lack thereof being perfectly understood. Nu- 
merical explosion models broadly match observations, even 
if they still cannot replace empirical models used to fit ob- 
served light curves. Cosmology is expecting to gain insight 
into possible distance biases from the steady improvement of 
numerical SN models. 

Several second-generation surveys have collected sam- 
ples of a few hundred spectroscopically confirmed and well- 
measured SNe la events up to z ~ 1 over the last decade. This 
world sample amounts in total to about 1000 events and the 



dark energy constraints already extracted from about one half 
of it ([242]) show that calibration and statistical uncertainties 
contribute similar amounts. Improving the photometric cali- 
bration, improving the survey design to make it less sensitive 
to calibration errors and increasing the statistics (in particular 
at the highest possible redshift) are the three main handles to 
improve over present constraints. Since spectroscopic iden- 
tification of several thousand events involves too much ob- 
serving time, broadband photometry should now also deliver 
identification of the events, which pushes for increased qual- 
ity measurements. Constraints from other dark energy probes 
will, hopefully, soon become available and crossing results 
will tighten and secure the obtained dark energy constraints. 

Future SN surveys are planned but only secured on ground- 
based facilities, although observing from space opens the pos- 
sibility of measuring precise distances to supernovae at z > 1, 
a key handle on dark energy constraints. Ground-based imag- 
ing surveys and space-based imaging and spectroscopy sur- 
veys target high-quality dark energy constraints. For the next 
decade, the LSST (ground) and Euclid (space) projects are 
aiming at surveying most of the extragalactic sky, and expect 
to see their first light around 2020. Currently, only LSST plans 
to constrain dark energy from supernovae distances. 
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